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ABSTRACT 


Results are reported of experiments on turbulent circular jets 
issuing into a cross flow, both for unheated and heated jets. 
Longitudinal and transverse distributions of velocity, pressure, 
temperature, and turbulence intensity are presented. 

The velocity distributions depend on the ratio of jet to 
cross flow momentum flux. The temperature distributions depend 
also on the density ratio. The jet structure is dominated by a 
vortex wake which forms behind the jet, which is evident both from 
the detail measurements and from smoke photographs of the flow. 

The experiments also indicate that the components of cross flow 
normal and parallel to the jet trajectory independently control 
the entrainment rate. 


ii 



ACKNOWL EDGMENTS 


The authors wish to acknowledge the inception of this work 
by the late Dr. H. K. Wiskind. They also acknowledge the financial 
support of the National Aeronautics and Space Administration through 
Grant NGR- 36-02 7-008. 


iii 



TABLE OF CONTENTS 


Page 

ABSTRACT ii 

ACKNOWLEDGMENTS iii 

TABLE OF CONTENTS iv 

LIST OF FIGURES vi 

CHAPTER I. INTRODUCTION 1 

CHAPTER II. GENERAL DISCUSSION 4 

1. Brief Review of Previous Work 4 

2. Coordinate System 7 

3. Dimensional Analysis 8 

CHAPTER III. EXPERIMENTAL EQUIPMENT 11 

1. Wind Tunnels 11 

2. Nozzle 12 

3. Air Supply System 14 

4. Probes 14 

5. Other Equipment 17 

CHAPTER IV. EXPERIMENTAL PROCEDURE 18 

1. General 18 

2 . Unheated .let Experiments 19 

3. Heated Jet Experiments 28 

CHAPTER V. RESULTS OF MEASUREMENTS 31 

1. Unheated Jet Experiments 31 

2. Heated Jet Experiments 35 

CHAPTER VI. DISCUSSION OF RESULTS 40 

1. Jet Centerline 40 

2. Detail Structure of Jet 46 

3. Entrainment Analysis 61 

iv 



CHAPTER VII. FLOW VISUALIZATION EXPERIMENT 73 

1. Introduction 73 

2. Equipment 74 

3. Experimental Procedure 75 

4. Results 76 

5. Discussion of Results 76 

CHAPTER VIII. SUMMARY AND CONCLUSIONS 81 

REFERENCES 84 

FIGURES 86 

TABLE 158 

APPENDIX A Pressure Force Acting on Jet in Cross Flow 159 

APPENDIX B List of Symbols 162 


v 



LIST OF FIGURES 


Page 


Figure 

i. 

Coordinate System 

86 

Figure 

2. 

Sketch of Wind Tunnel 

87 

Figure 

3. 

Diagram of Test Apparatus 

88 

Figure 

4. 

Structure of Nozzle 

89 

Figure 

5. 

Initial Velocity Profiles of Jet 

90 

Figure 

6 . 

Behavior of Free Jet from 1/4 in. nozzle 

91 

Figure 

7. 

Air Supply System 

92 

Figure 

8. 

Measuring Probes 

93 

Figure 

9. 

Constant-Current Hot Wire Set 

94 

Figure 

10. 

Examples of Hot Wire Calibration Curve 

95 

Figure 

11. 

Diagram of Flow Direction Coordinates 

96 

Figure 

12. 

Normal and Parallel Velocity Components 

97 

Figure 

13. 

Location of Jet Velocity Centerline 

98 

Figure 

14. 

Correlation of Jet Velocity Centerline 
Location 

99 

Figure 

15. 

Contours of Constant Velocity 

100 

Figure 

16. 

Contours of Constant Velocity 

101 

Figure 

17. 

Contours of Constant Velocity 

102 

Figure 

18. 

Contours of Constant Velocity 

103 

Figure 

19. 

Projection of Velocity Vectors 
on Normal Plane ( R =4) 

104 

Figure 

20. 

Projection of Velocity Vectors 
on Normal Plane ( =6) 

105 

Figure 

21. 

Projection of Velocity Vectors 
on Normal Plane ( R - 8) 

106 


vi 



Figure 

22. 

Rotational Velocity Field in Vortex 
Region ( ^ = 8) 

107 

Figure 

23. 

Rotational Velocity Field in Vortex 
Region ( £ =4) 

108 

Figure 

24. 

Rotational Velocity Field in Vortex 
Region ( =6) 

109 

Figure 

25. 

Rotational Velocity Field in Vortex 
Region ( =8) 

110 

Figure 

26. 

Contours of Constant Normal Velocity 
( ^ = 4) 

111 

Figure 

27. 

Contours of Constant Normal Velocity 

( R = 6) 

112 

Figure 

28. 

Contours of Constant Normal Velocity 

( K - 8) 

113 

Figure 

29. 

Mass Flux in Jet 

114 

Figure 

30. 

Distribution of 2 -Momentum Flux 

115 

Figure 

31. 

Distribution of Velocity along Jet 
Velocity Centerline 

116 

Figure 

32. 

Total Pressure Distribution 

117 

Figure 

33. 

Static Pressure Distribution 

118 

Figure 

34. 

Distribution of Turbulence Intensity 
( K = 4) 

119 

Figure 

35. 

Distribution of Turbulence Intensity 

< R - 8) 

120 

Figure 

36. 

Energy Spectra of Turbulence in Jets 

121 

Figure 

37. 

Centerline Location of Heated and 
Unheated Jets 

122 

Figure 

38. 

Location of Temperature Centerline 
( ^ = 4) 

123 

Figure 

39. 

Location of Temperature Centerline 

( R - 8) 

124 

Figure 

40. 

Contours of Constant Temperature 

125 


vii 



Figure 41. Contours of Constant Temperature 126 

Figure 42. Contours of Constant Temperature 127 

Figure 43. Contours of Constant Temperature 128 

Figure 44. Contours of Constant Temperature 129 

Figure 45. Contours of Constant Temperature 130 

Figure 46. Contours of Constant Temperature 131 

Figure 47. Contours of Constant Temperature 132 

Figure 48. Distribution of Maximum and Centerline 

Temperature 133 

Figure 49. Distribution of Centerline Temperature 134 

Figure 50. Temperature Distribution in Far 

Downstream Region 135 

Figure 51. Diffusion of Temperature in Jets 136 

Figure 52. Diffusion of Temperature in Jets 137 

Figure 53. Diffusion of Temperature in Far 

Downstream Region 138 

Figure 54. Comparison of Jet Centerline Location 139 

Figure 55. Correlation of Temperature Centerline 

Location 140 

Figure 56. Potential Core Coordinates 141 

Figure 57. Potential Core Length of Jet in Cross 

Flow 142 

Figure 58. Temperature Distribution in the Plane 

of Symmetry 143 

Figure 59. Distribution of Entrainment Velocity 144 

Figure 60. Mass-Momentum Flux Diagram 145 

Figure 61. Mass and Momentum Flux in Cross Flow 

Direction 146 

Figure 62. Temperature Centerline Location in 

Far Downstream Region 147 


viii 



Figure 

63. 

Smoke Generator 

148 

Figure 

64. 

Photographic Arrangement 

149 

Figure 

65 . 

General View of Jet in Cross Flow 

150 

Figure 

66. 

Close-up View 

150 

Figure 

67. 

Cross-Sectional View of Jet 

151 

Figure 

68. 

Cross-Sectional View of Jet 

152 

Figure 

69. 

Cross-Sectional View of Jet 

153 

Figure 

70. 

Cross-Sectional View of Jet 

154 

Figure 

71. 

Vortex Motion in Laminar Jet 

155 

Figure 

72. 

Coordinate System for Flow around 
Parabolic Cylinder 

156 

Figure 

73. 

<2 -Momentum Flux Distribution 

157 


ix 



I . INTRODUCTION 


In recent years considerable attention has been given to the 
problem of jets issuing into a subsonic cross flow. This is mainly 
because it is the basic configuration in many practical situations. 
In general the point of main interest is different in each case. 
Main fields of application and their associated problems can be 
summarized in the following list. 

1. V/STOL aircraft use air jets to obtain lift for take-off. 
During transition flight these jets are bent rearward, which 
changes the air flow around the aircraft body, causing a loss 
in lift and instability due to pitching moment. Choosing the 
number of jets and their locations are the most important 
problems in order to get maximum performance of the aircraft. 

2. A hot stream of gas in a pipe or in a duct can be cooled 
simply and effectively by injecting jets of cool gas into it. 
In order to get a desired temperature distribution at a given 
section downstream of the jets, it is necessary to study the 
mixing of two flows. 

3. A plume leaving a chimney with substantial upward momentum 
is often bent by a cross wind (bent-over plume) . In this air 
pollution problem it is important to calculate the ground 
level concentration for a given situation. It is customary 

to apply a suitable diffusion formula from a certain point at 
which the upward momentum of the plume is negligible. The 
analysis of a jet in a cross flow can be applied to find the 
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flow structure near this point as well as its location. 

4. A river with some contaminant materials, flowing into the 
ocean or the lake with cross current, presents another problem 
of diffusion. In this case the main interest is in finding 
the polluted region and the degree of pollution (concentra- 
tion distribution). 

The present study was motivated by the cooling problem. The 
temperature of hot combustor gas for a gas turbine is generally 
too high for the material used, and it needs to be cooled. Gener- 
ally, cooling is done by injecting cool air through holes in the 
surrounding wall. Although multiple jets are used for effective 
cooling, a single jet case is treated here as a first step. The 
results obtained give basic information not only for the combustor 
cooling problem but also for many other applications. 

Quite a few papers have appeared in the last decade concerning 
the experimental and the analytical studies of a jet in a cross 
flow. They studied mainly the over-all behavior of the jet such as 
the jet centerline location. In Chapter 2 these previous works are 
summarized together with some basic discussions. Compared with 
such general aspect of the jet, its detail structure has not been 
thoroughly investigated. Therefore, it was felt that instead of 
going into the cooling problem directly, it was necessary to 
investigate first the detail aspect of the jet injected into a 
cross flow of equal temperature, which is essential in order to 
understand the mixing process. 



3 


To show the jet structure in detail, distributions of velocity, 
pressure and turbulence intensity were measured at several sections 
of the jet. The data revealed a complex structure of the jet. The 
structure was found to be strongly influenced by a pair of counter- 
rotating eddies which were formed behind the jet. Thus the rota- 
tional velocity field was also investigated in detail. 

After this unheated jet experiment, the mixing of two flows 
with different temperature was studied. To do this, the jet was 
heated above the free stream temperature instead of cooled as in 
the cooling problem, because experimentally it was easier to get 
large enough temperature difference by heating. Although buoyancy 
force acts in the opposite directions in heated and cooled jets, it 
does not change the general flow characteristics significantly. The 
jet was heated up to 400°F. Lateral distributions of temperature 
at several jet sections and the longitudinal distribution along the 
centerline were measured in detail to study the mixing of two flows. 
All of the experimental results are presented in Chapter 5. The 
discussion based on these data is given in Chapter 6. 

In addition to these experiments, a flow visualization experi- 
ment using kerosene smoke \?as conducted. The pictures presented in 
Chapter 7 show clearly the cross-sectional shape of the jet in a 


cross flow. 



II. GENERAL DISCUSSION 
1. Brief Review of Previous Work 

Since the present flow is the basic configuration in many prac- 
tical situations, numerous papers have been published. The problems 
they treated can be classified mainly into three groups. The first 
group concerns the determination of the jet centerline location 
either experimentally or analytically. The second group treats the 
problems relating to V/STOL aircraft. The third group investigates 
theoretically the rolling-up of a jet due to a cross flow. 

Several forms of the jet centerline equation have been sug- 
gested. Some of these are empirical, and others are obtained by 
semi-empirical analyses. The latter usually contain several 
numerical constants which have to be determined from the experi- 
mental data. Some of the early works have been mentioned in Ref. 1. 
The recent work has been reviewed in Ref. 2 and 3. Some of these 
works will be mentioned elsewhere in this report. Several jet 
centerline equations reviewed in Ref. 2 show quite a scatter (see 
also Fig. 54 of this report). Carlson et al. (Ref. 3) tried to 
correlate some of the previously published data taking into account 
the different wind tunnel sizes relative to the jet sizes in the 
different experiments. 

To analyze the jet deflection due to the cross flow, two 
important factors have to be considered. One is the aero-dynamical 
force acting on the jet, and the other is the entrainment of the non- 
turbulent main stream into the jet. Most of the analytical models 
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consider only the former factor (Refs. 1, 4, 5). As the entrainment 
idea developed in recent years, it was realized that the effect of 
the entrained momentum on the jet deflection is more important 
except for the region close to the nozzle where the aerodynamical 
force is just as significant (Ref. 6). A few investigators obtained 
the jet centerline equation by considering the entrainment effect 
only (Ref. 7). 

In applications to V/STOL aircraft, the most important factor 
is the pressure distribution on the plate surface. Experimental 
data for the surface pressure distribution were given in several 
papers (e.g. Refs. 8, 9, 10). Theoretical analysis of the problem 
has also been presented by some investigators. The calculation 
made by Wu et al. (Ref. 10) was based on the simple two-dimensional 
potential flow model, and it gave reasonably satisfactory agreement 
with the measured wall pressure distribution. 

Recently theoretical analyses of the jet rolling-up process 
have been presented by several authors (e.g. Ref. 11). Their 
analysis is an extension of the early work by Chang (Ref. 12). 

Chang considered the time-dependent deformation of the cross-section 
of the hypothetical circular fluid column in a cross flow. In this 
analysis the original jet is replaced by a number of vortex fila- 
ments distributed along the jet surface; by following the loca- 
tion of these filaments as time goes on, instantaneous cross- 
sectional shape can be obtained. Since such analysis shows only 
the potential flow aspect of the interaction of two flows, the 
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result cannot be compared with the data which are taken usually in 
a turbulent flow. 

Perhaps most directly related to the present work are the 
experiments of Keffer and Baines (Ref. 6), Callaghan and Ruggeri 
(Rep. 13, 14, 15), and Jordinson (Ref. 16). Keffer and Baines made 
measurements in the plane of symmetry in the upstream region 
( for unheated jets. Their experiments tend to show 

the importance of the momentum ratio in determining the jet trajec- 
tory. They also studied entrainment, making an analytical model 
which is a simple extension of the free jet model. This is satis- 
factory because they were only concerned with the small I" region. 
Callaghan and Ruggeri made measurements using a heated jet. Their 
measurements also were only in the plane of symmetry, and their jet 
was rather closely confined because of their small wind tunnel width 
(3.1 to 7.8 nozzle diameters). Jordinson measured total pressure 
distributions for the full lateral extent of the jet, not only in 
the plane of symmetry. Recent developments in most of the problem 
discussed above are reported in Ref. 17. 

To provide an adequate understanding of the behavior of a jet 
in a cross flow, the present work examines the detailed lateral 
distributions of velocity, temperature and turbulence intensity, 
both in the upstream and downstream regions of the jet. Longitudinal 
distributions of the above quantities, ideas about entrainment, and 
the role of the vortex motions are deduced from the extensive 
lateral measurements. 
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2. Coordinate System 

Figure 1 shows the coordinate system adopted in the present 
work. The Cartesian coordinate system ( X , ^ , Z ) with the origin 
at the center of the nozzle exit is convenient for practical appli- 
cations. The 3C -axis is the direction of the cross flow, the 
Z -axis is the direction of jet injection, and the ^ -axis is 
perpendicular to both. On the other hand for the purpose of analy- 
tical study, it is more convenient to use the coordinate which fol- 
lows the deflected path of the jet. For this, the natural 
coordinate system based on the mean streamlines is most suitable. 
However, it is very difficult to determine this system experi- 
mentally. The coordinate system ( ^ ^ ) adopted here is some- 

what similar to the natural one. ^ is measured along a jet 
velocity centerline, which is one of the mean streamlines in the jet. 
The plane normal to this streamline at each point of ^ is assumed 
to represent the cross-section of the jet. Therefore the £ -axis 
(jet centerline) can be determined experimentally by connecting 
points of maximum velocity at each section in the plane of 
symmetry (JC^-plane with ^ =0). On each normal plane the ^ —axis 
is taken to be in the direction of the ^ -axis from the jet center 
(thus y - ), and the ^ -axis is perpendicular to the £ -axis. 
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3. Dimensional Analysis 

Because of the complicated configuration of the problem a 
purely theoretical analysis is very difficult, and it is beyond 
the scope of the present study. However, it is very useful for 
the experiment to find important dimensionless numbers by simple 
dimensional analysis. 

Consider a heated (or cooled) jet issuing normally to a uni- 
form flow • As the flow develops the jet becomes turbulent. Some 
mean quantity (Q. at some point downstream of the nozzle can be 
considered to be a function of the following quantities: 

G = (X, V, Z, P y 6 Uj, Po, Pj, f-7S, /?„ /£ Cf) 

where subscripts 0 and J denote the properties of the cross flow 
and the jet, respectively. The free stream and the jet are con- 
sidered to be the same fluids (no subscript for , Qp ) . 

In this analysis, only the case in which both the initial jet and 
the cross flow are laminar as in the present experiment is discussed. 
The effect of the surrounding wall is assumed to be negligible. 

Then important dimensionless groups formed from quantities listed 
above are expected to be as follows. 

S / x & Z ViD U,P Cf/A 
1 P ' P ' P ' P ' P > £ , 

Vi" A J4. jS_ ) 

P$o\ f-T.j J r. j o, ' p a J 
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The length scale for the cross flow has been chosen as p because 
the cross flow produces wake behind the jet whose charactertic 
length is P near the nozzle, the region where the wake begins 
to appear. 

Assume 



The last two assumptions lead to the conclusion in a free turbulent 
flow that the viscosity is not important, that is the Reynolds 
stresses are much larger than the mean viscous stresses. Then we 
obtain 

JL =/.* Hill ^ ) 

Q t + ( » ' p ' r. ' u. J r*' MPiTj-T*\\ <» 

\z 

number which expresses the ratio of the inertia to the buoyancy 
force in the initial jet. 

When 77-77 ( E - P. , since in this analysis 

the two fluids are assumed to be th® same kind), equation (1) 
is simplified to 



The quantity 




h 77 


called the Froude 


( 2 ) 
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Usually n - r so that in this case the only important 

quantity is the ratio of two velocities. 



III. EXPERIMENTAL EQUIPMENT 

The main equipment consists of a wind tunnel, an air supply 
system, a nozzle and measuring probes. An overall sketch of the 
wind tunnel is shown in Fig. 2, and the details of the setup are 
shown in Fig. 3. When designing the setup, care was taken to insure 
the following: 1) uniform, low-turbulent free stream. 

2) low-turbulent jet with very uniform and constant exit velocity 
distribution. 3) sensitive probes for velocity and temperature 
measurements. 4) accurate probe positioning device. 

1. Wind Tunnels 

For preliminary measurements a small wind tunnel with 
4" x 4" test section was used. It is an open-circuit tunnel with 
a centrifugal fan driven by a 1 HP DC motor. Using this wind tun- 
nel, the effects of test section walls and the tunnel wall bound- 
ary layers on the structure of the jet were investigated. This 
tunnel was also used for flow visualization experiments. 

The main experiment was conducted in a larger wind tunnel 
which has a 28" x 28" test section and a length of 15 ’6'' (Fig. 2). 
This tunnel is an open-circuit , suction type tunnel with an axial 
fan driven by a 53 HP DC motor. This wind tunnel had been used in 
the Jet Propulsion Lab. at Case Institute of Technology, and was 
reconstructed for the present experiment. It is a low-speed, 
low- turbulence wind tunnel, and thus very suitable for this experi- 
ment. Fine mesh screens are installed at the intake and serve to 
make larger eddies into smaller ones which decay faster. Due to 
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these screens and a high contraction ratio (10.4:1), the turbulence 
level at the test section is very low (0.4% at 20 ^ 30 ft /sec, the 
velocity range in which present data are taken) , and the velocity 
profile is very flat. The walls of the test section are made of 
Plexiglas plates. One side wall has a small door for easy 
access to the equipment inside the test section. 

The jet was injected at a position three feet downstream from 
the beginning of the test section. The boundary layer thickness 
at this position was 1-1/2 in. In order to minimize the effect 
of the tunnel wall boundary layer, the flat plate, made of Plexi- 
glas 17 in. wide and 24 in. long, was placed at a position 2-1/2 
inches above the bottom wall, and the nozzle is set flush xdLth 
this plate. The plate has a wedge-shape nose which is connected 
smoothly to the flat part. At the injection point, which is 2 
inches from the nose, the boundary layer thickness is small com- 
pared with the diameter of the nozzle (less than 10%) so that its 
effect is expected to be negligible. 

For the heated jet experiments it is desirable to have con- 
stant free stream temperature. Since this is an open-circuit 
wind tunnel, temperature cannot be controlled and is subject to 
the outside temperature change. However, the tunnel was placed 
in a well air-conditioned room, so that the temperature change 
during any one run of the experiment was within + 0.2°F. 

2. Nozzle 

One of the important factors which control the structure of 
jet is the velocity profile at the exit of the nozzle. Generally 
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the velocity profile is not flat in practical cases; however in 
order to maintain well defined experimental conditions the initial 
velocity profile was made as flat as possible. For this purpose a 
specially designed nozzle was used in the present experiment. Its 
structure is shown in Fig. 4. Six turbulence-reducing screens 
of various mesh sizes are used, and the contraction part has 
a carefully designed contour. The contraction ratio is 16:1, 
and the exit diameter is 1/4". Fig. 5 shows velocity profiles 
at the exit for various jet speeds. Profiles are very flat over 
the whole diameter except for the thin boundary layer near the 
wall. Turbulence level at the exit is 0.3%. 

For the heated jet experiment the nozzle was wrapped with 
asbestos for thermal insulation. The temperature profile at 
the exit was also found to be very flat. 

To see the general behavior of the free jet from this nozzle, 
distributions of the mean velocity and the turbulence intensity 
along the jet centerline were measured, and are shown in Fig. 6. 

It can be seen that the length of the potential core (in which 
the flow is laminar, and the velocity remains constant) is four 
times the nozzle diameter. The location of the jet virtual origin 
(the origin of the self-similar region) is found to be very close 
to the nozzle exit (0.4 times the nozzle diameter upstream of 
the nozzle) . Beyond ^/p = 3 , the centerline velocity de- 

creases hyperbolically with the distance from this virtual origin, 
and the velocity profile is expected to be similar. 
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3. Air Supply System 

Air for the jet was supplied by a 100 HP, two-stage compressor 
which can supply 402 CFM of air. Since air pressure of the com- 
pressor varies in the range between 80 and 100 PSI, it was 
necessary to design a suitable air control system in order to 
get constant velocity at the exit of the nozzle. The schematic 
of the system is shown in Fig. 7. Pressure fluctuation is 
reduced first by the pressure regulator to the order of 
+ 1/4 in. of water. This fluctuation is further reduced by the 
flowrate regulator so that the fluctuation of the flow rate is 
within + 0.5% of the set value. The flowrate was measured by a 
Fisher and Porter rotameter. 

For the hot jet experiment, air was heated by a heating 
tape wrapped around a steel pipe of 1-1/2 in. diameter. Using 
this heating tape it was possible to make jet temperature at the 
nozzle exit up to 400°F at the maximum flowrate. 

4. Probes 

Hot Wire and Anemometer 

The hot wire probe is made from brass tubing of 1/4 in. dia- 
meter and 26 in. long with wire supporting tips made of bronze 
needles (Fig. 8). The diameter of the hot wire is 0.00015 in., 
and it is made from 90% platinum-10% rhodium Wallas ton wire. The 
wire is etched in a thin column of nitric acid. The etched por- 
tion is about 0.045 in. long. The end of the hot wire probe is 
connected by a Cannon connector to the wire which leads to the 



15 


input of the hot wire anemometer set, A constant-current hot wire 
set, Shapiro and Edwards Model 50, was used in the experiment. 

The set consists of an amplifier, resistance bridge, potentiometer, 
mean square output meter and square wave generator, as shown in 
Fig. 9. The potentiometer allows one to measure the wire current 
up to one hundredth of a milliamp, which makes accurate velocity 
measurement possible. The time constant of the hot wire was 
about 450 microseconds, so that for the measurements of turbulence 
intensity and energy spectrum, the compensation circuit was 
applied to compensate forthe thermal inertia of the wire. 

Pressure Probes 

Total pressure was measured by a Kiel total head probe which 
is insensitive to pitch and yaw angles up to about 40 degrees. 
Where space is limited near the nozzle, a total head tube made 
from a hypodermic needle of 0.028 in. diameter was used (Fig. 8). 
Static pressure was measured by a Pitot static tube made from a 
hypodermic needle. It has two static pressure holes of 1/64 in. 
diameter at the positon 7 times tube diameter downstream from a 
rounded nose. 

In a turbulent flow both total and static pressure measure- 
ments are subject to errors due to the lateral turbulent motion. 
Errors are not small in the present experiment, and this problem 
is discussed in the next chapter. 

Temperature Probe 

Temperature was measured by an Iron-Cons tan tan thermo-couple 
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supported by a brass tube of 1/4 in. dia. and 26 in. long (Fig. 8). 
The diameter of the tip contact is 1/16 in. 

Y aw Angle Meter 

For the measurement of the mean flow direction, a small sensi- 
tive yaw angle meter was used. This is made of two hypodermic 
needles of 0.028 in. diameter with bevelled open ends. The flow 
angle can be determined by rotating the meter until pressures 
measured by the two tubes become equal. To measure three- 
dimensional flow direction, usually four tubes are used. But this 
increases the size of the probe and the rotating device becomes 
more complicated. In the present experiment the flow angles 
were measured in two planes to get three-dimensional flow direction. 

In a turbulent flow the measurement is affected by the Reynolds 
stress, which is discussed in the next chapter. 

Probe Positioning Apparatus 

In the experiment most data were obtained by point by point 
measurements using various probes. Thus it was necessary to use 
an accurate probe positioning apparatus. The one used was a 
three-dimensional manipulator manufactured by Brinkmann Instrument. 
It can traverse 4, 3, and 3 in. in each direction, and can rotate 
the probe in the vertical plane. In order not to disturb the 
flow, the positioning device was put outside of the test section 


as shown in Fig. 2. 
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5. Other Equipment 
Micromanome ter 

Pressure was measured by a Meriam micrometer Model 34 FB2 TM 
reading in thousandths of an inch of water pressure. 
Microvolt-ammeter 

The thermo-couple was connected to a microvolt-ammeter 
(Keithley Instruments, Type 150A) . Since this is a sensitive 
instrument, it is very difficult to measure the mean temperature 
in a highly turbulent flow, because the instrument follows the 
fluctuations. To damp out large fluctuations R-C circuit was 
connected to the meter. 

Wave Analyzer 

Energy spectre^ of turbulence were measured by a General Radio 
wave analyzer Model 1900A. This has a frequency range of 
0 to 55,000 cps and three bandwidths of 3, 10 and 50 cps . The 
input was taken from the amplifier of the hot wire set. Accuracy 
of the measurement is + 3% of the indicated value up to 50 KC. 



III. EXPERIMENTAL PROCEDURE 


1 . General 


The present experiment consists of two parts. The first one 
is the investigation of an unheated jet in a cross flow, in 
which distributions of velocity, pressure and turbulence intensity 
are measured. The next one is for a heated jet, in which the 
temperature distribution is measured in detail. 

In the unheated jet experiment, the most important variable 
is the velocity ratio /Uo . In the heated jet case we hgve 
one more variable Va . In many occasions these two 


variables are coupled and appear as a single parameter Uo 

which is called the momentum ratio. Since most of the previous 


works treated the unheated jet cases ? the term 

( i 

effective velocity ratio R which is the square root of the 
momentum ratio is used in this report for the purpose of making 
it easy to compare the data in both experiments. 

The value of R can be chosen arbitrarily, but to cover the 
practically important range, three cases R 4, 6, 8 were 
investigated inboth heated and unheated jet experiments. These 
nominal values of R correspond to using the free stream speed 
upstream of the plate as the cross flow speed. Since the free 
stream speed above the plate is a few percent higher than that 
ahead of the plate, the actual values of R appropriate to the 
experiment are slightly lower than the nominal values. The nominal 
and actual values of R , along with the corresponding free stream 
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and jet speeds for the unheated jet experiment, are as follows. 


R (nominal) 

4 

6 

8 

R. (actual) 

3.91 

5.80 

7.73 

do (ft/sec) 

30.7 

25.8 

20.7 

dj (ft/sec) 

120 

150 

160 


For Simplicity , the nominal values will be quoted hereafter. The 
above combinations of speed were chosen so that differences in 
both jet and cross flow speeds for three cases were small. 

For the heated jet experiment, only jet speeds were varied 
according to the jet temperature (density) keeping cross flow 
speeds constant. Experimental conditions of this case are sum- 
marized in Table I. Settings of both heating current and flow rate 
meter were determined by trial and error at .first to get desired 
values of P* Ul . For high temperature jets it took about 

c t d 

two hours to achieve an equilibrium temperature. Care was taken 
to have constant jet velocity temperature and cross flow speed 
throughout any one run of the experiment. Detail procedures for 
two cases (unheated and heated) are explained below. 

2. Unheated Jet Experiments 
Jet Centerline 

First, the location of jet centerlines was measured for 
/\ = 4, 6, 8. In addition, the experiment was done also for 
^ - 10. These profiles are important not only for practical 
reasons but also for defining the coordinate system. Experimentally 
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the centerline profile can be obtained by connecting points of 
maximum velocity in the plane of symmetry at each section of the 
jet. As shown in the next chapter the velocity is not necessarily 
maximum at the jet center of each section. This method of deter- 
mining the jet centerline is inevitably trial and error, because 
the normal section of the jet cannot be specified until the center- 
line profile is known. In the present experiment normal planes 
to the centerline were determined by a step by step method in which, 
starting from the horizontal plane near the nozzle (jet was injected 
vertically upward), each normal plane was approximated by extra- 
polating the previously determined centerline. 

The point of maximum velocity at each section was found both 
by the hot wire and by the total head tube (Kiel probe) . Both 
results agreed very well. Since the velocity along the centerline 
decays much faster than the one in a free jet it is very difficult 
to determine the jet centerline beyond 25. 

Velocity Distribution 

For each momentum ratio, the mean velocity distributions were 
measured at several sections of the jet by the hot wire. Since 
the mean velocity changes its direction mainly in the Xi^-plane 
(the plane parallel to the cross flow) , the hot wire was set 
normal to this plane. It is known that a hot wire measures the 
velocity component normal to it. The mean flow direction 
measurement (which is explained later in the chapter) shows that 
the change of the flow direction in the JCW -plane is relatively 
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small except in the small region where the rotational velocity 
is large. Therefore the velocity measured in the manner described 
above is expected to be close to the magnitude of the actual velo- 
city (the difference is about 3% for 15 degrees change in the 

-plane) . 

The velocity was measured point by point traversing the hot 
wire according to the coordinate defined above. In order to see 
the detail distribution, it was necessary to measure at 60 'v* 100 
points for each section depending on the actual size of the jet. 
This was a very cumbersome procedure, but since the actual size 
of the jet was small, it was impossible to insert a rake without 
disturbing the flow very much. An automatic traversing device 
in an inclined plane was not easy to design, and was given up 
for the present experiment. 

In the present procedure it took about two hours to finish 
one section. Therefore the important problem was the 'drifting' 
of the hot wire calibration curve on which the velocity calculation 
was based. This drifting was made small by warming the wire 
for a rather long time (about 1.5 hours) just before the data were 
taken. Also, in order to minimize the error during the experiment 
the drifting of the wire was checked in a certain interval by 
measuring the free stream velocity, and the calibration curve 
was corrected accordingly. The overheat ratio of the wire was 
set at 0.3. Care was taken to have a very straight calibration 
curve in the operating xange. Fig. 10 shows the examples of 
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the hot wire calibration curve. 

The velocity measured by the hot wire is not exactly equal 
to the actual velocity in a highly turbulent flow. The relation 
between the measured and the true values is given by the equation 
(Ref. 18), 

-/. / ) 

U ( / _ T U 1 u 1 -' 

where (A is the longitudinal velocity fluctuation and IT is the 
later/al fluctuation. This correction was not done in the ex- 
periment, because the turbulence level is rather low in the jet^ 
as shown in the next chapter. 

Me an Flow Direction 

Using a small yaw meter the mean flow directions were 
measured at points where the mean velocities were known. As has 
been explained in the previous chapter, this yaw meter was made 
of two hypodermic needles so that it can determine the flow 
direction only in one plane. Therefore in order to measure 
the three-dimensional flow direction as in the present case, it 
is necessary to measure in two planes which are perpendicular to 
each other. In the present experiment flow directions were 
measured in the XZ -plane and the -plane. To determine the 

mean velocity vector at each point, calculations are necessary 
using the magnitude of the velocity measured by the hot wire 
and flow directions in two perpendicular planes. Referring to 
Fig. 11, their relations are given by 
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where notations should be clear from the figure. When (3 is 
small, U k Um as explained in the previous section. 

The purposes of the flow direction measurement are to cal- 
culate the mass and the momentum fluxes across each section of 
the jet, and also to obtain the distribution of the rotational 
velocity (the velocity component parallel to the normal plane). 
Thus it is necessary to know velocity components normal and 
parallel to the jet cross section. They are given by (Fig. 12) 


u n =• U* (<*'- 0 ) 

Uf> = (4 [a 4 V + £0rp - C&d (o(— ojJ 
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where Q is the angle of inclination of the jet centerline 
(or the value of <X at the center) , and (f is the angle between 
Up and the parallel component of the cross flow velocity. 
Although the flow in the jet has been assumed to be parallel to 
the jet centerline at each section, in the region away from the 
center the flow becomes parallel to the free stream. Consequently, 
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in order to see the rotational velocity field in the vortex region 
(which is of major interest here) one has to project the velocity 
on the plane normal to the path of the vortex core. This can be 
done simply by replacing Q in Eq. (4) by the angle determined 
from the path of the vortex core. 

It can be shown (Ref. 19) that the behavior of the yaw meter 
is expressed by the equation 



where and are pressures measured by two tubes of the 

meter, Q is the angle between the mean flow direction and the 
yaw meter, and {/IT is the Reynolds stress. The above equation 
shows that even if 0 » 0, p, and f? are not necessarily 
equal in a turbulent flow with shear. Since the turbulent jet 
is a typical shear flow, this yaw meter gives an error depending 
on the value of the Reynolds stress at the measuring point. In 
a free jet the maximum value of - is about 0.02, so that 

the maximum error in @ is about one degree. In the present 
experiments the local mean velocity is much larger than the one 
in the free jet because of the convective velocity (J 0 . Thus 
the error in is expected to be much smaller than the above 
estimate. 

Since this procedure of mean flow direction measurement is 
more painstaking than the velocity measurement, in order to save 
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labor it was decided to measure the flow direction only in the 
half of the normal plane, assuming that the flow is symmetric with 
respect to the center-plane. The equipment was carefully set up 
so that the flow was almost symmetric, as shown in the next 
chapter. 

Since to calculate the mean velocity vector it is necessary 
to measure three quantities (speed and two angles) at one point, 
care was taken to place the probe at exactly the same point for 
each measurement. 

Pressure Distribution 

At present there is no reliable way to measure static 
pressure in a highly turbulent flow. An ordinary Pitot static 
probe suffers appreciable error due to lateral velocity fluctua- 
tions. The correction of the measured value to obtain the true 
value of the static pressure is complicated in the sense that it 
depends on the lateral turbulence intensity and also on the 
ratio of the average eddy size of turbulence to the diameter of 
the Pitot tube. In the present experiment another problem is 
the alignment of the probe, since the mean flow direction changes 
very much at different points. Therefore the static pressure 
was not measured, but was calculated from the values of the total 
pressure and the speed. A Kiel total head probe was used for the 
total pressure measurement. This probe is particularly convenient 
when the exact flow direction is unknown. Similar to the Pitot 
static probe, the measured values should be corrected in a 
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turbulent flow. The correction formula is given by 





/? 
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where both and u/- are laterial velocity fluctuations. Since 
!/^ and were not measured, the measured values of /^. were 
not corrected. But the error is expected to be small 
because the turbulence level is rather low in the flow due to 
the convective velocity 
Turbulence Intensity Distribution 

The distribution of the turbulence intensity in the normal 
plane was measured by the hot wire. Since the wire was set normal 
to the X.Z -plane as in the velocity measurement, the turbulence 
intensity measured in this way gives the fluctuation level in the 
direction of . However it is expected to be very close to 

the intensity in the direction of the mean velocity vector 

The constant-current hot wire system which was used in the 
present experiment is not appropriate for measurement in a highly 
turbulent flow, because of the non-linearity of the hot wire 
response. It is supposed to be accurate in a flow with the 
turbulence intensity less than 10%. As shown in the next chapter 
the maximum intensity observed was 22% in the present case, but 
in the most of the region the intensity was less than or slightly 
larger than 10%. Error is 10 ^ 15% if the turbulence intensity 


is about 20%. 
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The wire overheat ratio was 0.3 as in the velocity measure- 
ment, and the x^ire response was compensated for thermal inertia 
in order not to lose the energy of the high frequency fluctuation. 
Energy Spectrum 

The turbulence energy spectra were measured at a few points 
in one section of the jet, and they were compared with the spectrum 
at the corresponding position in a free jet. 

The output of the hot wire was first amplified and compensated, 
then sent to the wave analyzer. The mean square output voltage 
of the wave analyzer at frequency is given by 




where is the filter bandwidth and ) is the 
energy/cps at j* = . The actual spectrum E(f)o f the 
velocity fluctuation is calculated from the relation 


E (f) — 


<$(f) 
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where Q is the transfer function of the hot wire set, and 
is the sensitivity of the hot wire. 

Since 


/ E (+) df 
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then 

s (/j _ fi(-k) _ Sffi. (7) 

^ ~ei uT (j-) 

This is a convenient form of presenting the energy spectrum. The 
filter bandwidth A f of 3 cps was used in the experiment. 

3. Heated Jet Experiments 
Jet Velocity Centerline 

As in the unheated jet measurements, the jet velocity center- 
line location was measured by the hot wire and the total head probe. 
When the jet temperature is very high, the effect of density changes 
cannot be neglected. Since the hot wire responds to the 'mass 
velocity' (J , the mean electric current through the wire (for 
a constant current set) corresponds to the value of f (J . Separat- 
ing velocity and density into mean and fluctuating parts, we have 

Ju — T U -+ T 7 * . 

Therefore the point of maximum U does not necessarily correspond 
to the point of the maximum velocity. Neglecting the correlation 
between density and velocity fluctuation which is very small usually, 
it is still necessary to know the mean density (or temperature) distri- 
bution in order to find the point of maximum (J . It was found that 
this correction was very small because the density change was smaller 
than the velocity near the center. In this respect the measurement 
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by the total head probe is better, because it can detect the variation 
of n A (the static pressure is almost uniform in the jet). How- 
ever, since the sensitivity of the manometer was limited, the total 
head probe was used mainly in the upstream region. 

Temperature Centerline 

The point of the maximum temperature in one section does not 
coincide with the maximum velocity point, so that one must define 
separately the centerline for the temperature distribution in a 
manner similar to the jet velocity centerline. Temperature was 
measured by the IC thermo-couple. The temperature centerline 
location was obtained by connecting points of maximum temperature 
in the plane of symmetry at each jet cross secjxon which is already 
specified by the jet velocity centerline. 

Temperature Distribution 

The temperature distributions were measured at several 
sections of the jet for - 4 and 8. The data were taken at 
60 ^ 100 points in one section, and they were connected smoothly 
to get uniform contours of constant temperature. The thermo- 
couple was connected to the micro-volt meter to measure the small 
temperature difference. The fluctuation of the reading was made 
small by the use of the additional RC circuit. Temperature 
distribution at each section was measured up to the point where 
the temperature is equal to the mean of the maximum and the free 
stream temperature. Beyond this point the measurement becomes 
difficult due to the large fluctuation. The temperature at the 
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nozzle exit was checked and kept constant throughout any one 
run of the experiment. 



V. RESULTS OF MEASUREMENTS 


1. Unheated Jet Experiments 

Fig. 13 shows the jet velocity centerline location for the 
effective velocity ratios ^ = 4, 6, 8 and 10. These profiles 
can be correlated by a single curve 


JL 
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as shown in Fig. 14. The constant in the right hand side represents 
the effect of the potential core. The power for Ua /u • was 

d 

found by trial and error. The power for and other numerical 

constants were determined by drawing smooth fitting curves. 

The error of the measurement increases as we go downstream, because 
the velocity distribution becomes very flat near the center. The 
estimated error in the centerline location near the most downstream 
section measured was + 3%. 

Using the coordinate system determined from the jet velocity 
centerline obtained above, the distributions of the velocity 
at several jet sections for ^ = 4, 6, 8 were measured, and 
results are shown in Figs. 15 through 18. Strictly speaking, 
these are the distributions of Um instead of (J . Although 
( J can be calculated from (J^ according to Eq. (3), the difference 
is small, and it is not important in order to observe the jet 
structure. The results are presented in the form of contours 
of the non-dimensionalized velocity difference U which is defined as 
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U — Uo 

^ Uo 


Although the use of (/ is just a matter of convenience, at far 
downstream sections where the jet velocity is almost parallel to 
Uo it is possible to classify the flow configuration as 




o< U ■£ / 

tc<o 


jet region 
wake region 


At upstream sections the directions of two velocities {J and 
are different so that it is not simple to distinguish between jet 
and wake. 

Since the velocity {J is composed of the axial velocity \J^ 
(normal to the jet cross section) and the velocity yj (parallel 
to the jet section), it is important to separate these two compon- 
ents. After the mean flow directions were measured in two planes 
at each point, Up and cT were calculated according to Eq. (4). 
The distributions of Up are shown in Figs. 19, 20, 21. The data 
were taken in one half of the plane, because the flow is almost 
symmetric with respect to the centerplane. Fig. 22 shows the 
rotational velocity field in the developing region (of the 
rotational motion) for =8. Figs. 23, 24, 25 show the rota- 
tional velocity field at downstream sections for = 4, 6, 8. 

In these figures the velocities \j and are non-dimensionalized 
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by , which appears to be the appropriate scale, because the 

rotational motion is induced by the cross flow. 

One of the important features in the study of a turbulent 
jet is the analysis of the mass flux increase which is directly 
related to the entrainment of the surrounding non-turbulent fluid 
by the turbulent flow. In order to calculate the mass flux through 
the normal plane at each section, it is necessary to know the dis- 
tribution of the normal velocity . Uy) can be calculated 

from Urn according to Eq. (4). Then similar to 7T > the dimen- 
sionless quantity (J ^ for the normal velocity is defined as 



Un (Jo a 


where Uo/) ~ cos Q is the normal component of the cross flow 

A/ 

velocity. In Figs. 26, 27 and 28 the contours of constant Um 
are shown. In these figures only the jet region is considered 
for the purpose of the mass flux calculation. The use of the 
quantity is more convenient than the previously defined , 

A'' 

because does not include the rotational velocity whereas [X 
is related to both a^ial and rotational velocities. The calcula- 
tion of \J^ requires the measurement of the flow direction only 
in one plane, as can be seen from Eq. (4), which simplifies 
measurements and calculations compared with those of Up . This 
made it possible to measure the distributions of Uy) at a 
larger number of sections than were possible for Up , and 

measurements of U at several sections are necessary to study 

n 
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the mass flux change with distance. In order to study the mass 
flux increase in a turbulent region from the mean velocity dis- 
tribution, it is necessary to define the T edge f of the turbulent 
region. Following Keffer and Baines (Ref. 6), the jet boundary 
is defined as the surface where the normal velocity excess 
( Ur, - U.„ ) is 10% of the maximum value at a given section; 
that is, the turbulent region is where o.K K/iil* Then the 
mass flux can be calculated by integrating \J h over this region. 
The momentum flux can be calculated in a similar way. Fig. 29 
shows the result of the mass flux calculation. For the entrainment 
analysis in Chapt. 6, the momentum fluxes in the X and Z 
directions are calculated by decomposing (Jn into both directions. 
Fig. 30 shows the Z -momentum flux distribution. 

The velocity distribution along the jet centerline is shown 
in Fig. 31. The ordinate was chosen as the inverse of the non- 
dimensionalized (by (Jj ) velocity (J^ in order to see clearly 
the hyperbolical decrease of (Jc near the nozzle. For comparison, 
free jet data are shown in the same figure. The data were taken 
in the free from the same 1/4 in. nozzle with the initial velocity 
Uj = 120 ft/sec. In each case the virtual origin (of the flow 
near the nozzle in cases of the jet in the cross flow) can be 
found by extrapolation of the straight part as done in the figure. 

The distributions of the total pressure were measured at a 
few sections just to compare with Jordinson f s data (Ref. 16). 

One of the present results is shown in Fig. 32. Notice that the 
data were taken not on the plane normal to the jet centerline, 



but on the plane normal to the cross flow ( -plane), similar 

to Jordinson. The static pressure distribution at the same posi- 
tion is shown in Fig. 33. The values of the static pressure were 
calculated fo^rm the values of the total pressure and the velocity 
measured by the hot wire. At several points the static 
pressure was measured by the Pitot static probe with its axis 
aligned in the mean flow direction at each point. Both results 
for the static pressure agreed within the experimental error. 

Figs. 34 and 35 show the distributions of the turbulence 
intensity relative to the maximum velocity at the section 

for R> = 4 and 8. The hot wire gives the quantity f(?7 \j 
(intensity relative to the local velocity) , however the distribu- 
tion of LA 2 " itself is more interesting because it shows the 
distribution of the kinetic energy of turbulence. 

To see the more detail structure of turbulence, the energy 
spectra were measured both at the center of the jet and at the 
point of maximum velocity (center of the vortex motion) at the 
section Vp = 23 for =8. The results are shown in Fig. 36. 
The energy spectrum of turbulence in the free jet at the corres- 
ponding position is also shown in the same figure for comparison. 
The figure shows the slightly different structure of turbulence 
in two jets. 

2 . Heated Jet Experiments 

To see the similarity of the jet structure with respect to 
the momentum ratio, the jet velocity centerline locations were 



36 


measured in the heated jets for ^ = 4 and 8. The result is 
presented in Fig. 37. The figure shows that the velocity center- 
line location of the heated jet coincides with that of the 
unheated jet as long as the momentum ratio is kept constant. The 
initial jet temperature above the free stream temperature was 
75 °F and 320°F. The velocity ratio and the density ratio for 
each case are listed in Table 1. 

Figs. 38 and 39 show the profiles of the temperature center- 
line (locus of the maximum temperature in the plane of symmetry) 
of heated jets in comparison with jet velocity centerline profiles 
for /\ = 4 and 8. From these figures it is clear that the tem- 
perature centerline location cannot be correlated by the momentum 
ratio alone. 

Similarly to measurements of velocity distributions, tempera- 
ture distributions were measured at several sections both in the 
mildly heated jet and in the high temperature jet. The results 
are presented in Figs. 40 through 47. They show the contours of 
the dimensionless coefficient 9 which is defined as 



The location of the maximum temperature at a given section is 
not necessarily in the plane of symmetry as is also true for the 
velocity distribution. Unlike velocity, temperature is a scalar 
quantity so that only one dimensionless temperature has to be 
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defined. The above choice of Q is convenient, as is seen to 
vary from 0 to 1 at any one section. 

One important use of the measurements of temperature distri- 
butions is that it is possible to investigate the jet structure at 
the far downstream position where it is no longer possible to 
measure the velocity distribution. In the present experiment the 
data were taken at several positions up to = 70 for the high 
temperature jet ( — /© = 320°F). For this jet the effect of 

buoyancy cannot be neglected. It is possible to see the effect 
on the jet structure by comparing temperature distributions of 
the highly heated jet with those of the mildly heated jet. This 
is discussed in the next chapter. For the large value of £ , 
which is measured along the jet velocity centerline, ^ was 
determined by extrapolating the centerline which was neasured 
up to about ^/p = 25. The normal plane for large ^ was 
taken to be the vertical plane normal to the cross flow -plane) , 

since the angle of inclination of the centerline was very small 
in that region. 

Fig. 48 shows the distribution of the maximum and the center- 
line temperature along ^ in jets with the initial temperature 
77-7: = 320°F for — 4 and 8. For ^ = 4, the maximum 

and the centerline temperature are identical. 

The figure also shows the temperature distribution in the 
free jet with "77 ~ To = 330°F. It is known that in a free jet 

7 

the temperature excess ( 7£<yr To ) decreases inversely with 
the distance from the virtual origin just as velocity. For the jet 
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in the cross flow the same behavior as the free jet is expected 
for small ^ where the jet still retains the characters of the 
free jet. In Fig. 49 the inverse of the dimensionless temperature 
excess is plotted as the ordinate so that the hyperbolic change 
can be seen clearly. Also in this way it is easy to find the 
virtual origin. The higher temperature curves in Fig. 49 are 
simply replots of the data of Fig. 48. 

Fig. 50 is a plot similar to Fig. 49, but for the downstream 
region. Here one finds that ( 7m* - To ) </->*). 

This means that the distribution of the temperature along ^ can 
be expressed by the relation 



where ^ represents the location of the virtual origin of the 
temperature distribution for large ^ , which is different from 

the one found for the distribution for small ^ 

From distributions of temperature at several sections of 
the jet, one can see the lateral spread of temperature with dis- 
tance. To show this quantitatively, the half-value radius lfy z 
was introduced. is defined as the effective radius of 

the area ( S'/ z ) occupied by the contour of & * 0.5 where the 
temperature excess ( T-To ) is 50% of its maximum value 
( 'Ttob'X ~To ) at a given section, namely 
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is analogous to the half-value radius defined for the 
circular jet. When the temperature distribution is similar at 
each section as in a free jet. the use of as a characteris- 

tic length scale in the lateral direction is very convenient. In 
the present case the temperature distribution is almost similar 
so that the above defined is expected to represent the 

lateral spread of temperature in the jet. Fig. 51 shows the 
distributions of fk with the distance ^ in jets with 

- J~o = 320°F for R - 4 and 8. Free jet data are also 
shown. In fig. 52 the distribution of in jets with dif- 

ferent initial temperature are shown. The heat spreads linearly 
near the nozzle as in the free jet. However, in the far downstream 
region it seems that it spreads with a certain power of ^ 

This power was found in the same way as the temperature distribu- 
tion along £ . Fig. 53 shows that in the far downstream 
region the temperature spreads as 



where the value of is not exactly equal to but close to 

the value given in Eq. (9). 



VI. DISCUSSION OF RESULTS 


1. Jet Centerline 

According to the dimensional analysis in Chapter II, there 
are two non-dimensional numbers, namely I Uo and Jj / f* 0 , 

which are important in the problem of a jet in a cross flow, when 
the buoyancy force is negligible. From a physical consideration 
it is known that the deflection of the jet is predominantly 
determined by one dimensionless quantity, y* I ft (the 

momentum ratio). 

In order to check the importance of the momentum ratio experi- 
mentally, it is necessary to conduct an experiment using a heated 
(cooled) jet or using two different types of fluid. However, if 
two densities differ very much, the buoyancy force becomes important 
and complicates the analysis. 

Eq . (8) in Chapter V which was obtained by correlating the 
data taken in the present unheated jet experiment shows the 
importance of the velocity ratio. From the heated jet experiment, 
however, it is clear that the momentum ratio is indeed the 
significant factor to determine the jet deflection (Fig, 37). 
Although the effect of buoyancy is not negligible for the high 
temperature jet, the centerline location can be correlated by 
the momentum ratio alone, at least in the range of temperature 
used in the experiment. One then expects that the deflection 
angle of the jet is determined mainly by the ratio of the jet 
momentum flux in the Z -direction to the momentum flux in the 
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?C -direction. 

When the buoyancy force is not negligible, the jet centerline 
tends to shift upward (if K fc )• But this tendency is com - 
pensated by the increase of the momentum flux in the 7C -direction 
resulting from the increased rate of entrainment due to the 
decrease of the jet density (this is discussed more in a later 
section). Therefore, it is difficult to distinguish the effect 
of buoyancy from the study of the jet centerline location. 

Eq. (8) can be written in the form (even when /*, ^ f* 0 ) 
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where z yp = 0,^-2 R represents the point from where the 

deflection of the jet starts, and is often called the effective 
source of the jet (not to be confused with the virtual origin of 


the jet). For 6, this point nearly coincides with the end 


of the potential core. Since for /?^4, the potential core 
itself deflects, the two points do not coincide. Thus, if R 
is increased, Z& should approach to a value equal to the length 


of the free jet potential core, so that the above relation be- 
tween j£ 0 and R is true only in a certain range of R . The 
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general case is discussed in the next section. For practical pur- 
poses consideration of the effective source is not convenient. 

For this reason most empirical formulas of the jet centerline 
set the origin at 2 = 0. If the jet deflection is assumed to 
start at 2 ~ 0, the present data give the following equation: 


( 11 ) 


Comparison of Eq. (8) (or Eq . (11)) with other published 
data is presented in Fig. 54 , These data were taken by Jordinson 
(Ref. 16), Keffer and Baines (Ref. 6), Margason (Ref. 2), Storms 
(Ref. 20) and Ivanov (Ref. 1). The empirical equation of the 
centerline obtained from those data is given in the form 

| = A ^ ( jrf 

where /\ , ^ are numerical constants. The value of " 2. " 

given by the above investigators ranges from 0.67 to 1.0, and 
the value of " " from 0.33 to 0.39. This wide scatter of 

data is attributable to the differences in experimental conditions. 
The major factors seem to be the initial conditions of jet and 
the wind tunnel size. The jet initial conditions (velocity 
profile, pressure, turbulence level) are not given in most of 
the papers. The sizes of the wind tunnels used by the above 
investigators are summarized Ref. 2. According to the investiga- 
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tions by Carlson et al. (Ref. 3), the wind tunnel size in the pre- 
sent experiment as well as in others is large enough so that its 
effect on the jet path is negligible. From the analysis of the 
free jet it is known that the jet structure becomes independent 
of the initial condition after some distance downstream from the 
nozzle. However, the centerline location of the jet in a cross 
flow at a given section is affected by the upstream condition. 

So far there is no systematic study about the effects of the jet 
initial conditions on the jet path. Since these conditions in 
the present experiment are close to the standard ones (flat 
velocity profile, low turbulence level), the centerline location 
is about the middle of other data as can be seen in Fig. 54. 

Due to the complexity of the problem, there is no purely 
theoretical calculation of the centerline location available. 
Several semi-empirical formulas have ben derived. Vizel and 
Mostinski (Ref. 4) calculated the deflection of the centerline 
from the ratio of the jet momentum flux in the Z -direction 
(assumed constant) to an estimate of the pressure force acting 
in the X. -direction. McAllister (Ref. 7) used a similar method, 
but instead of the estimated pressure force he used an estimate 
of the momentum flux in the X -direction which is brought into 
the jet by the entrained fluid. In addition to special assumptions 
contained inthe estimates, both analyses finally require that a 
pair of numerical constants be determined from the experiments, 
which of course makes it easier to fit the derived formulas to 
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the data. 

As for the temperature centerline, its location is different 
from that of the jet velocity centerline. As can be seen in 
Figs. 38 and 39, unlike the jet velocity centerline location the 
momentum ratio alone cannot correlate the temperature centerline 
location. For a given momentum ratio, increasing jet temperature 
causes less penetration of the temperature centerline. In order 
to find the important factors which control the temperature field, 
it was tried to correlate the temperature centerline location 
of heated jets with various combinations of velocity and density. 
The results are presented in Fig. 55. The temperature centerline 
location is found to be represented by a single curve 



This combination of density and velocity ratios does not mean 
much physically. It is necessary to find the dimensionless 
quantities which are more important in the problem. Since the 
temperature field is strongly related to the velocity field, the 
momentum ratio is expected to be one important factor to determine 
the temperature centerline location. As mentioned above the 
momentum ratio alone cannot correlate the data. To find another 
quantity, one has to understand why the temperature centerline 
shifts with increasing jet temperature while the jet centerline 
is fixed (for a given momentum ratio) . As is explained in 
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detail in the next section, when the jet temperature is increased, 
the strength and the size of the vortex motion are increased. 

Thus, more heat is transferred into the vortex region, and more 
concentration of heat around the vortex center results. The 
temperature centerline does not mean much physically but it 
merely represents the over-all penetration of heat into the cross 
flow. Therefore its location is strongly affected by the location 
of the vortex center, especially when much heat is concentrated in 
the vortex region. Thus, as the jet temperature is increased for 
a given momentum ratio, the vortex center moves away from the jet 
center (due to the increased size of the vortex motion) , and con- 
sequently the temperature centerline shifts downward. Such 
increased activity of the vortex motion seems to be attributed 
to the increased rate of entrainment. Ricou and Spalding (Ref. 21) 
found that in a free jet the rate of entrainment is related to 
the density ratio £ //. . Thus it is expected that increasing 

jet temperature (decreasing jet density) causes more entrainment, 
which activates the vortex motion. Besides the jet density, 
the velocity (J ' and the kinematic viscosity are changed 

< 7 </ 

when the jet temperature is changed while keeping the same momentum 

ratio. As has been discussed in the dimensional analysis 

(Chapter II), the flow structure is assumed to be independent 

of the viscosity since the Reynolds number is very large. As 

for U; , it is increased when 77 is increased. But its 
<7 d 

change is rather small (inversely proportional to the square 
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root of the density change), and its effect on the vortex motion 
seems to be less important compared with the density effect. 

The buoyancy force also increases the vortex motion. However, 
this seems to be a secondary effect in the present experiments, 
because the temperature difference is already small at the 
position where the buoyancy force becomes important to the vortex 
motion. Consequently Eq. (12) is better written in the form 



Unlike the jet velocity centerline, the temperature centerline 
starts deflecting from a point very close to the nozzle, so that 
it was unnecessary to find an effective source of the coordinate 

corresponding to 2? 0 for the jot centerline. 

2. Detail Structure of Jet 
Velocity Field 

Right at the exit of the nozzle the jet is considered to be 
potential flow, because of the low turbulence intensity and very 
thin boundary layer. Due to the very steep velocity gradient 
at the edge of the jet, the flow is very unstable, which causes 
the development of the free shear layer (called mixing layer). 

The mixing layer thickness increases with distance, and the 
potential flow region (called potential core) disappears after 
some distance from the nozzle. In the potential core the velocity 
remains essentially constant. Except for the small momentum ratio 
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( R £ 4) , the potential core is deflected very little by the 
cross flow. For a free jet this core length depends on the jet 
Reynolds number, boundary layer thickness at the exit, and 
turbulence level. For a jet in a cross flow it has been found 
from the experiments that the core length decreases with decreasing 
momentum ratio. The reason for this is that the mixing layer 
spreads faster when there is cross flow, because of the tangential 
stress induced by the cross flow. The relation between the core 


length and the velocity ratio (for simplicity, the case 
% - /» is considered) can be obtained as follows. First, 

consider the free jet. The relation between the spreading angle 
of the mixing layer £ 0 and the core length £ 0 is given by 
(see Fig. 56) 



A jet in a cross flow is subjected to the additional shear 
stress which acts tangentially to the jet cross-section. The 
increase of the spreading velocity of the mixing layer (entrainment 
velocity) due to this tangential stress is assumed to be proportion- 
al to the characteristic tangential velocity scale, which is 
chosen to be the cross flow velocity. Then the increase of the 
spreading angle is proportional to Uo / (Jj ■> since the boundary 

<r 

between the mixing layer and the potential core is moving with 
a velocity nearly equal to (Jj . Denoting the additional 
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spreading angle by £ ' which is assumed to be much smaller 
than , we obtain 

d 

where Q is a constant. Therefore the core length is given by 


/ / 

p (u) 

Z u< 

d 

This relation is compared with the data taken by Pratte and Baines 
(Ref. 22) in Fig. 57. In their experiment ^°/p - 6.4 and Q, - 0.36. 
As can be seen in the figure Eq. (14) agrees very well with the 


data. Kronauer (Ref. 23) obtained exactly the same relation by 
correlating his experimental data. He gave = 6.0 and 

C = 0.67. 


After this initial region, the jet becomes fully turbulent, 
and the jet deflection commences. In this region of the flow 
the important features are the distortion of the jet cross-section 
and the formation of the vortex wake behind the jet. The distor- 
tion is mainly due to the pressure distribution around the jet and 
the entrainment. The jet spreads to the side of lower pressure. 
Since this spread-out region has smaller axial momentum than the 
region near the center, it is deflected more by the pressure drag 
and the entrained momentum in the cross flow direction. Consequently 
the jet cross-section h^s the shape of a crescent. This deforma- 
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tion is also influenced by the vortex motion in the wake when 
it becomes strong. It is well known that behind a cylinder placed 
normally to a uniform stream a pair of counter-rotating 1 standing 
eddies 1 are formed in the Reynolds number range 6 K 
If the Reynolds number is increased more, eddies shed from the 
cylinder alternately, forming the familiar Karman vortex 
street in the wake region. Very similar to this situation, 
a pair of eddies are formed behind the jet soon after the jet 
leaves the nozzle. These eddies do not shed, however, because 
of their strong interaction with the jet; instead they acquire 
axial momentum from the jet, and move along the jet path. 

These eddies are somewhat similar to a 1 vortex-pair T . After the 
formation the development of the vortex motion is controlled by 
such factors as the normal component of the cross flow-velocity, 
the shape of the jet cross-section (which acts to block the cross 
flow) and the jet initial velocity (not including the case 
J5+fo ). For a given velocity ratio (for a given jet 
path and a given jet cross-sectional shape), increasing the cross 
flow velocity causes faster development of the vortex motion, but 
this is compensated by the increased jet velocity (to keep ^ 
constant) which gives less time for the development. Therefore, 
the development of the vortex motion is expected to be controlled 
by the velocity ratio. As the velocity ratio is increased, the 
jet experiences a larger normal component velocity of the cross 
flow at a given distance from the nozzle, because the jet is 
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deflected less. Consequently the strength of the vortex is 
larger for larger ^ (Figs. 23, 24, 25). Initially the vortex 
strength increases with distance. However, viscosity acts to 
retard this motion so that the vortex strength reaches its maxi- 
mum at a certain distance from the nozzle, and starts decaying 
thereafter. 

Fig. 15 shows that for ^ = 4 the jet cross-section is not 
much altered by the vortex motion, and that this vortex wake 
strongly interacts with the ordinary wake (which is formed by 
the cross flow separation from the jet in the region near the 
nozzle) since the height of the jet centerline is low. For larger 
R , the jet is swept around into the wake region due to the 
strong vortex motion (Figs. 16, 17, 18). 

After the initial region the jet still retains the general 
characteristics of a free jet for some distance, as evidenced 
by the hyperbolical decay of the velocity along the jet centerline 
(Fig. 31). But, the vortex strength grows rapidly, and its in- 
fluence on the jet structure increases. This region may be 
called the transition region. The jet centerline has very large 
curvature in this region. 

In the far downstream region, the jet centerline becomes 
nearly parallel to the cross flow, and the velocity excess in the 
jet (the velocity defect in the wake) is very small. The small 
velocity excess makes it difficult to study the jet structure 
from velocity measurements. However, it is possible to investigate 
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this region by measuring temperature distributions in the heated 
jet, or by using flow visualization techniques. 

Near the end of the transition region, the vortex motion 
nearly ceases developing, but it decays much slower than the axial 
velocity excess so that the vortex motion dominates the flow field 
in the far downstream region. 

Although the vortex motion is weak for ^ = 4, it is strong 
enough to retain the crescent-shaped jet cross-section for quite 
a long distance, and then as the vortex motion decays, this 
shape is destroyed gradually by diffusion. For larger momentum 
ratios , more axial momentum is convected by the rotational velocity 
toward the vortex region from the region around the jet center. 
Consequently the point of maximum velocity shifts from the jet 
center to the vortex center after some distance from the nozzle. 

By the time the development of the vortex motion is almost finished, 
most of the axial momentum is in the vortex region. As the vortex 
activity decays , this concentrated momentum in the vortex region 
diffuses in the turbulent region due to turbulent transport, 
and the distribution becomes more uniform. 

Temperature Field 

In the initial region the temperature at the center remains 
constant. The core length of the temperature is shorter than that 
of the velocity because the thermal boundary layer is a little 
thicker than the velocity boundary layer at the nozzle exit. At 
the end of initial region, the temperature distribution has very 
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flat contours. 

In the transition region, at first the centerline temperature 
decays inversely with the distance measured from the virtual 
origin, and then decays more slowly (Fig. 49). Due to the 
diffusion of temperature into the wake region, the centerline 
temperature decays faster than the temperature in a free jet. In 
spite of the deformation of the jet cross-section, the temperature 
spreads linearly at first as in the free jet, and then more 
slowly. This situation is similar to the jet in a co-axial uniform 
stream. When the jet temperature is increased (the jet density 
is decreased) , the centerline temperature decays faster because 
of the increased entrainment of the outside cold fluid (Fig. 49). 
The effect of the jet temperature on the jet spreading is not so 
simple. Increasing the jet temperature means increasing the 
jet velocity to keep the momentum flux f? (J; constant 

& V 

c fit 4 is kept constant in these experiments at any one 
momentum ratio). Therefore although the lateral spreading rate 
is increased, so is the axial velocity. This is the reason why 
the effect of the jet temperature change on the jet spreading is 
not so apparent in Fig, 52 as on the temperature decay, shown in 


Fig. 49. 

Contours of constant temperature have more rounded shape 
than those of constant velocity, since temperature diffuses more 
readily into the wake region than the momentum. This difference 
between the turbulent transport of heat and momentum can be seen 
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clearly for ^ = 4. As the momentum ratio increases, more heat 

concentrates around the vortex center. The relation between the 
degree of such concentration and the vortex activity is quite 
apparent for /(* = 8 when temperature distributions at several 
sections are compared. Up to about ^/p = 25, the heat concentra- 
tion in the vortex region continues to increase (Figs. 40, 41). 
According to the rotational velocity measurement (Figs. 22, 25), 
at Vp =23 the vortex activity has already passed its peak (the 
rotational velocity is larger at = 12.5 than at Vp = 23). 

After the peak up to around Vp = 25, the rotational velocity 
is still large enough to be able to resist the diffusive action 
of turbulence and keeps convecting heat into the vortex region. 

Then as the vortex motion decays, the heat is diffused, and the 
temperature distribution becomes more uniform in the turbulent 
region (Fig. 42). 

For a given momentum ratio, as the jet temperature increases 
more heat is convected into the vortex region (compare Figs. 41 
and 44), which indicates the increased vortex activity. This 
is because of the increased entrainment, as has been explained 
in the section on the jet centerline. 

In the far downstream region constant temperature contours 
become circular for K = 4, but later they become slightly 
oval-shaped, which shows the existence of the vortex motion at 
very far downstream positions (Figs. 46, 47). For ^ = 8 as 

explained above the balance between the turbulent transport and 
the convection by the vortex motion changes the temperature distri- 
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but ion. 

It is interesting to note that in the far downstream region 
the maximum temperature excess ( ^ ) changes with dis- 

tance as^Eq. (9)j 

7fi)\X To 


S-T. 


oC 




-0.77 


Also, the temperature profile spreads as^Eq. (10)^ 
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The configuration of the flow in the far downstream region is 
qualitatively similar to the free jet in a coaxial uniform stream. 
For a jet in a parallel stream it can be shown from the similarity 
analysis (Ref. 18) that the asymptotic variations of ( ~7to*>X. - To ) 
and are given by 



(15) 


Compared with this, the temperature in the jet in a cross flow 
decreases slightly faster, and also the jet spreads a little 
faster. This slight increase of the spreading rate appears to be 
due to the vortex motion which increases the entrainment rate. 
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Despite the gradual change of temperature distribution around the 
center of the turbulent region as explained above, the main jet 
structure remains almost similar as evidenced by the close similarity 
to the jet in a parallel stream. 

In the far downstream region, the lateral spread of tempera- 
ture in the plane of symmetry is predominantly determined by the 
axial convection and the lateral diffusion. As the rotational 
motion decays, the lateral convection by the rotational velocity 
becomes negligible compared with the axial convection. In the 
direction normal to the plane of symmetry ( -direction) the 
temperature is relatively uniform, so that diffusion in that 
direction is assumed to be negligible. Therefore, in analogy to 
the jet in a parallel stream, the temperature distribution in the 
plane of symmetry is expected to remain similar in the far down- 
stream region. 

Mathematically, the balance between the axial convection and 
the lateral diffusion is expressed as 

where P T is the eddy diffusivity in the % -direction, and £ 
is measured from the temperature center. Using the results 
obtained above, one can assume 

J - j r = #?) ~ T-Z = f(T) r* 

7c 4 
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where anc ^ ■n are constants. Substitution of these relations 
into Eq. (16) gives 

O' '$r+*r) 07) 


Putting ^ = 2 #2. “ 1* Eq- (17) becomes 

i 


(18) 


If one can assume £ = const , — 0 as in cases of a free 

jet and a jet in a parallel flow, Eq. (18) becomes 


Cf"+ W Uo Z f / -+ w Uo f — O 


(19) 


Using boundary conditions, 


/=/ 
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-f 0 
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the solution of Eq. (19) is 


i- TiH 

Tc-Z 


- y-ft V. 


( 20 ) 
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The values of n and C a re expected to be a function of the 
momentum ratio and the density ratio. However, Eq. (20) shows 
that by using a suitable lateral length scale, it is possible 
to obtain a universal temperature profile, and that the profile 
is very much like the Gaussian error curve. 

The experimental data presented in Fig. 58 show that this 
is indeed the case. The data were taken in the far downstream 
region ( i/p >/S ) for = 4 and 8 with ( Tj -T ) 

ranging from 50°F to 320°F. ^ ^ and £ are the values 

of (measured from the temperature center) at which is 

equal to 1/2 ( 7 1 * To ). The figure shows that in the 

positive region temperature falls slightly faster than the 
Gaussian curve near the outer edge. In the negative ^ region 
the lateral temperature distribution is almost Gaussian. This 
slower decrease in the negative £ region is due to the limited 
space between the jet and the bottom plate. In any case, the 
temperature distribution seems to be represented by a single 
curve. Practically this result is important because the tempera- 
ture distribution in the plane of symmetry gives the common 
basis for comparing the over-all mixing of a hot (cool) jet and 
a cross flow. 

Turbulence 

Measurements of velocity and temperature distributions give 
very useful information about the structure of the jet. However, 
it is important to understand that both distributions are very 
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much influenced by the turbulence structure of the jet. It is 
known that in a free turbulent flow the transport of quantities 
such as momentum, temperature and turbulence energy are related 
to the double structure of the turbulence (large eddies and small 
energy-containing eddies). The large eddies whose scales are 
comparable to the width of the flow, distort the bounding surface 
betyeen the turbulent and the non-turbulent flow, and thus control 
the entrainment process. These large eddies contain smaller 
scale eddies. Most of the turbulence energy is distributed among 
these small eddies (energy-containing eddies). The turbulent 
transport by the large eddy motion is considered to be convective, 
whereas the transport by small eddies is the gradient-type diffu- 
sion. The turbulent transport of a certain quantity can be 
influenced predominantly by either one of those features of 
turbulence, or by both. In the present experiments the vortex 
motion in the jet appears to be large eddy motion. It can be 
seen by comparing distributions of velocity and temperature 
(especially for ) that the transport of temperature is 

more strongly influenced by the vortex motion (or large scale 
eddies) than the transport of velocity. The transport by large 
eddies tends to make the distribution of any quantity uniform in 
the turbulent region, as evidenced by the large uniform tempera- 
ture region inthe middle for ^ =4 (Fig. 47). For large 

momentum ratios the velocity distribution is affected by the 
vortex motion, but this must be considered differently from the 
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large eddy activity. When the strength of the vortex is very large, 
it is no longer large eddy motion; rather it is a part of the mean 
flow field. 

Figs 34 and 35 show that the distribution of turbulence energy 
is very much like that of temperature, which means that the trans- 
port of turbulence energy is also strongly influenced by the con- 
vective motion of the vortex wake. It is known that the turbulence 
energy is nearly constant inside the turbulent region except near 
the bounding surface. Therefore the measured time-mean value of 
the turbulence intensity at any point depends mainly on this con- 
stant value and also on the intermittency factory Y ; namely we 
have 

? = r f: 

where is the measured turbulence kinetic energy per unit 
volume, and is its value in the turbulent region. The tur- 

bulence is assumed to be almost isotropic in the present experiments 
so that H 3 - . Since is known to be almost constant 

in the turbulent region, the distribution of (Figs. 34, 35) 

nearly represents the distribution of the intermittency factor. 
Generally Y is unity well inside the turbulent region.. In the 
present case Y is expected to be almost unity at the point when 
is maximum, and this point is considered to be the center of 
the turbulent region (which includes both jet and wake) as opposed 

to the previously defined jet center which has not as much physical 
meaning . 
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It is important to measure the turbulence energy distribution 
among eddies of different size for studying the structure of the 
isotropic turbulence, because the theory can predict some form 
of the energy spectrum in a certain range of the wave number which 
can be tested by the experiment. The energy spectrum is not 
measured very often in a turbulent shear flow, but it helps in 
understanding the detail structure of turbulence. Fig. 36 shows 
a compariosn of energy spectra measured in a free jet and in a 
jet in a cross flow for ^ =8. The spectrum of a free jet 
decreases according to - 5/3 power of frequency (or wave number) 
in a certain range of frequency. This region is called the inertial 
sub-range which is usually observed in isotropic turbulence with 
high turbulence Reynolds number, and this form of the spectrum 
can be explained by dimensional analysis for isotropic turbulence. 
Also, in a shear flow this - 5/3 power is observed in a rather 
wide range of frequency. This is mainly because eddies are 
usually isotropic in the higher range of frequency regardless of 
the anisotropicity of the energy-containing eddies (so called 
local-isotropy). The energy spectra of a jet in a cross flow 
measured at the center of the jet at the center of the vortex motion 
have slightly smaller power (- 1.75 power) than - 5/3 in the 
inertial sub-range. The spectral form of the energy containing 
eddies at the vortex center (which is very close to the center 
of the turbulent region as explained above) is very similar to 
the one at the center of the free jet, but the former decreases 
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faster in a higher wave number range, which indicates less dis- 
sipation of the turbulence energy. Compared to the spectrum at 
the vortex center, the one at the jet center has less energy 
in the energy-containing range, and more energy in the inertial 
sub-range. 

3. Entrainment Analysis 

An important feature of a free turbulent flow is the existence 
of a sharp boundary between turbulent and non-turbulent regions . 
Across such a boundary the random vorticity jumps from zero value 
in the non-turbulent region to non-zero value in the turbulent 
region (actually this vorticity change takes place in a very thin 
layer, the so-called viscous super-layer). Due to this sharp 
vorticity gradient, the vorticity diffuses into the non-turbulent 
region through the action of the viscosity. Once the vorticity 
is diffused, the strength of the vorticity associated with the 
velocity fluctuation is increased due to the stretching caused 
by the local mean rate of strain, and thus the originally non- 
turbulent region becomes fully turbulent . In this way the turbulent 
front moves forward into the non-turbulent region, and consequently 
the flow rate in the turbulent region increases with distance. 

This process is called entrainment . Since this entrainment pro- 
cess is very important in studying the structure of a free turbulent 
flow, it has been given much attention in recent years. It is 
known that the detail process of entrainment is influenced not 
only by the local flow conditions but also by the whole flow 
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structure. Thus its study involves all kinds of problems encoun- 
tered in the general study of turbulence structure. For this 
reason, the detailed entrainment process is not well understood. 
The recent stage of its study can be seen in Ref. 24. According 
to Townsend (Ref. 25) and others, it is possible to calculate 
the entrainment rates for the flows with structural similarity 
such as jets, wakes and turbulent boundary layers without re- 
ferring to the detailed process. Thus, the object of the present 
analysis is to study the entrainment rate for the jet in a cross 
flow assuming its similarity of the flow structure. 

In an axisymmetric flow the rate of increase of the mass 
flux in the turbulent region is given by 


Mi 


= 2.71 f 0 £ 


* df T * 

W 



( 21 ) 


where is the mass flux in the turbulent region (assumed to 

be non-dimensionalized) at a given section, 7j is the radius 
of the turbulent region, and a is the axial velocity at 
which the boundary is moving. Or using a simpler form 


dwtr 


2-Kfi yf u t 


( 22 ) 


where (Jg is the velocity at which the turbulent front is advanc- 
ing into the non- turbulent region, and it is called the entrainment 
velocity which represents the entrainment rate. The relation (22) 
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can be applied to the present problem as long as Yf is considered 

to be the effective radius such that = ( Siy where S' 7 - 

is the area occupied by the turbulent region. In the present 

experiment neither nor was measured, so that the 

velocity cannot be calculated directly from Eq. (22). The 

exact value of Ue. is not of interest, but its variation with 

distance is of interest in this analysis. Therefore, if it is 

v'""* 

possible to measure and r which are proportional to 

— - 

and T r ? respectively, at each section, the behavior of (J g can 
be studied. To do this, the jet boundary has been defined in 
Chapter 5 so as to measure WL and from the mean velocity 

distribution. 

In a circular jet the integration of the axial velocity over 
a cross section gives Mr itself, and in a region of similarity 
the lateral length scale such as the half-value radius for the 
velocity distribution represents the behavior of r T . For a 
circular jet in a coaxial uniform stream, the situation is not 
simple. The integration of the velocity over a given section 
is not convergent. The integration of the velocity excess con- 
verges, but its value does not mean much physically. In ths 
sense the situation is the same for wakes and turbulent boundary 
layers. Therefore it is necessary to define the boundary to get 
quantities related to Ift-p anc * *T 9 as ^ as been done in the 
present work. For the jet in a parallel stream the flow configura- 
tion near the nozzle is very similar to a free jet as long as the 
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jet initial velocity is very much larger than the free stream velo 
city. Thus, the measured 7% and (by defining the jet 

boundary) are expected to be proportional to W! r and r-F as 
in the free jet case. In the far downstream region where the velo 
city excess is very small compared with the free stream velocity, 
the velocity excess distribution is expected to have similarity 
structure. Since in this region the velocity is almost equal 
to Uo (free stream velocity) everywhere, 77L and are 

related by 

_m_ tc r**Uo _ , r* F 

' xr*^u„ ~ [ ~i f/ < 23 > 

inhere the ratio ' t Yj is expected to be constant because 

of similarity. In between the two regions there is no simple 

7 ^* it 

relation between Tfl and , Y and r T ■ 

In the present experiements , velocity data were taken in the 
regions corresponding to the intermediate and downstream regions 
of the jet in a parallel stream. Although the flow configuration 
of the present experiment is much more complicated than the jet 
in a parallel stream the velocity distributions show nearly 
similar patterns at each section. Thus, in analogy to the jet in 
a parallel stream the measured mass flux (shown in Fig. 29) 
is assumed to represent 70^ , and so is Yf . The entrainment 
velocity Ue calculated from measured M and f according to 
Eq. (22) is presented in Fig. 59. 
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In flows with similarity structure the entrainment velocity 
is known to be a constant fraction of the mean velocity variation 
at a given section, namely 

Ve = E ( U^ K - U 0 ) ( 24 ) 

where the constant of proportionality £ is called the entrainment 
constant. In these cases one can think of two fluids (turbulent 
and non- turbulent) flowing in the same direction. However, in 
the present case the two fluids are not parallel. Thus, it seems 
necessary to take into account the flow direction in the relation 
like Eq. (24). In fact, the value of £ calculated from measured 
and according to Eq. (24) is too large compared 

with the value for the free jet, and in spite of the nearly similar 
flow structure, £ changes very much with distance (it varies 
from 0.15 to 0.39 between = 7 and 23 for ^ =8, whereas 

in the free jet £ - 0.045). the important thing is that in 

the present problem the rate of strain field (which controls the 
entrainment process as explained above) is different from the 
case of two parallel streams. The rate of strain field should 
be decomposed into the parallel (to the jet path) component and 
the normal component. Notice that the length of the potential 
core of the jet in a cross flow can be obtained by the same 
consideration as explained in the previous chapter. Also, the 
normal component of the cross flow velocity is different from the 
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parallel part in the sense that it strongly influences the vortex 
motion. This vortex motion controls the bounding surface in a 
similar way to large eddies, and thereby plays a role in the 
entrainment process. Thus similar to Eq. (24) the entrainment 
velocity in the present problem is expected to be related to 
the velocity variations in two directions through the relation 

Ue = E/( Umw- U 0 OM 9 )i-E i (J 0 M*Q ' ( 25 ) 

The second term of the right hand side of the above equation 
represents the effect of the normal component of the free stream. 
Since it is known that the entrainment process is influenced by 
the whole flow structure, the two right hand side terms of 
Eq. (25) are not likely to be completely independent. This suggests 
that for a given velocity ratio (in this analysis only the 
case fj = f* 0 is considered for simplicity) the values of 
and jP ^ are nearly constant, assuming similarity of the flow- 
structure, but these values are different for the different velo- 
city ratios. 

Fig. 59 shows that Eq. (25) gives the variation of (Jq 
with very good accuracy using the constant values of and 

• The angle Q has been determined from the jet center- 
line profile. The values of ^ and ^ used in the figure 


are 
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R 

£, 


4 

0.070 

0.320 

6 J 

0.061 

0.242 

8 

0.067 

0.182 


These precise values may not be important, since Ue was not 

measured directly, nor calculated from Eq. (22) using VI j and 

h • The important fact is that constants fc: anc ^ /r 2 l 

exist. As can be seen in the velocity distributions for = 4, 

the wake region occupies the larger part of the turbulent region 
so that measured 'Tfl which is the mass flux in the jet (the 
region where (Jy ^ &cr)) is considered to be a small fraction of 
. However, for =8 the jet region occupies nearly 

the whole turbulent region. Thus, the values of and 

given above for =8 are probably close to the actual values. 

Eq. (22) can be considered as an over-all continuity equation 
across a jet cross-section. In addition to this we have momentum 
equations in two directions, namely in the axial direction 
( -direction) and in the normal direction ( -direction). 

Denoting the representative flow direction at a given cross- 
section by 0 (which is assumed to be equal to the angle of 
inclination of the jet centerline), and assuming that momentum 
changes only because of entrainment, the rate of change of 
momentum inthe axial direction is given by: 
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( 26 ) 


where M is the momentum flux across the jet section. Similarly 
the rate of change of momentum in the normal direction is given 
by: 


yu // * ■ „ (Jth 


(27) 


By solving the above equations, the following relations among 
n , M , and Q can be obtained: 


Do (M-Wr) = 


Mt J 

■&»$ 


M _ 

M, 

( ±Lf= 

[ Mr/ 


/ 


uj- 

H r 


!+ Th ( W-Wr) j 


(28) 


(29) 


(30) 


where and rir are the values of the mass and the momentum 

flux at the effective source of the jet. Eq. (30) which shows 
the relation between // and 'Tfl , is compared with the data in 
Fig. 60. As the momentum ratio increases, the curve approaches 
that of the free jet, which appears in the figure as a straight 
line M tMr - i. 
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In Eq. (28) the term ( ) represents the increase 

of the mass flux due to the entrainment, so that assuming the 
entrained fluid has the average velocity Uo , the term 
Uo ( *1 ~ Wf ) represents the increase of the momentum flux 
in the X -direction, namely 

Ml — Uo (7H-7r) r ) ' (31) 

In Fig. 61 this relation is checked using the experimental values 
of M x and n • The figure shows that for small Mr 

is substantially larger than the entrained momentum. The difference 
is due to the pressure force. In the far downstream region the 
effect of the pressure force seems to be negligible. 

In formulating Eqns. (26) and (27) the effect of the pressure 
field has been completely neglected. As a result of this assump- 
tion the momentum flux in the Z -direction is expected to remain 
constant without any force acting in that direction. As explained 
above, the momentum flux increase in the JC -direction is due 
to both the entrained momentum and the pressure force for small ^ 
Therefore, it is expected that the Z -momentum flux decreases 
especially in the region near the nozzle. The data show the 
gradual decrease of the -momentum flux with distance from the 
nozzle (Fig. 30). In that figure the absolute value of is 

not important, especially for = 4 and 6, because Mz 

is calculated only in the jet region neglecting the Z -momentum 
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in the wake region . Such decrease of Mz partly due to the 
definition of the jet region, but mainly due to the effect of the 
pressure. In the Appendix the pressure effect is estimated by 
considering the pressure distribution on the jet surface. 

Assuming that the 2 -momentum approaches a constant value 
after some distance from the nozzle where the pressure force is no 
longer important, the equation for the jet path is given from 
Eq. (28) as 


9 = 


dz 

dz 


Ht 

Uofa-VrJ . 


(32) 


Using the relation 




Eq . (32) becomes 


dz 




o 


(33) 


The boundary conditions for the above differential equation are 



where the effective source is assumed to be at the origin within 
the accuracy of this analysis (so that also /% = «<>). In 
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In order to solve Eq. (33) it is necessary to know 
which cannot be given at present by theory alone as may be clear 
from the expression of given in Eq. (22). Therefore, 

the analytical derivation of the centerline equation is always 
semi-empirical. Using experimental data for , Eq. (33) 

can be solved for small ^ and large ^ . Since there are several 
analyses available concerning the jet centerline location near the 
nozzle, Eq. (33) is not solved for small £ here. Instead, it is 
tried to solve Eq. (33) in the far downstream region using the 
information obtained from the temperature measurements. 

When £ is very large, we have ^ and Eq. (33) 

becomes 



(34) 


According to Eq. (10) the area of the turbulent region for large ^ 
0 72 

varies as XL Since the velocity is almost equal to (J Q 

everywhere in this region, the mass flux increases as Uo K °* 72 , 
that is 


w =- 


KP*U 0 (-£-) 


0.92 


a 


= K'P Uo (■*-) 


- o 2? 


(35) 
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where /C / is a constant, but its value is different for the different 
momentum ratios. Substituting Eq. (35) into Eq. (34), and using 
the boundary conditions, the solution is given by 


-2- — M. p 1 - / x ) a2f 
/> fc' * ( p y 


(36) 


where - % . The power of 0,28 is exactly equal to the 

value obtained by Pratte and Baines (Ref. 22) from the study of a 
flow-visualization experiment. In the present experiment, too, the 
profile of the temperature cneterline for large ^ has the power 
very close to 0.28 (Figs. 62). 



VII. FLOW VISUALIZATION EXPERIMENT 

1. Introduction 

In fluid mechanics as well as in other fields, visual observa- 
tion is the most convincing way of understanding what is really 
happening, especially when the situation is rather complicated. 
However, we are very often plagued by problems involved in making 
good flow visualization. The main problems are the selection of 
a proper tracer, the method of its injection into the flow, and 
obtaining good photographs. 

In the present experiment kerosene smoke was used as tracing 
gas, because it is easy to produce it continuously, and the smoke 
is clear to see. Although the smoke is not so toxic, it is 
necessary to have proper ventilation of the room. Another trouble 
was that kerosene smoke condensed everywhere in the pipe line, and 
this caused fluctuation of smoke pressure. This situation was 
improved by putting drains at proper positions. The smoke was 
injected into the same nozzle that was used in the main experiment. 

Although both side views and cross-sectional views of jet are 
interesting, cross-sectional views are particularly important to 
help elucidate the jet structure, and therefore cross-sectional 
views were taken at a number of downstream sections. Illuminating 
the cross-section of the jet by light from a thin slit, pictures 
of cross-sectional views were taken. The equipment used for this 
experiment, procedures and results are explained below. 
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2. Equipment 
Smoke generator 

The structure of the smoke generator is shown in Fig. 63. The 
smoke obtained from this generator is premixed with air so that 
the pressure and the flow rate of the smoke can be varied, which 
is convenient for the present experiment. Kerosene is heated 
above its boiling point in an evaporator which is made from a glass 
tube with a coiled heater. The kerosene vapor is then mixed with 
cool air, and becomes oversaturated. Thus small drops of kerosene 
are formed which can be seen as smoke. The flow rate of air can 
be adjusted and is measured by a flowmeter. Generally the 
increase of the flow rate after the air is mixed with kerosene 
vapor is negligible. The liquid kerosene level is kept constant 
in an evaporator so that the amount of smoke generated is controlled 
by the temperature of the heater. 

Wind Tunnel 

For convenience of taking pictures a small wind tunnel with 
4x4 in. test section was used. The test section is 22 in. long, 
and is made of Plexiglas plates. The turbulence level of the free 
stream was 0.3%. The jet was injected at a position 2 in. down- 
stream from the beginning of the test section. At this position 
the boundary layer thickness was about half of the nozzle diameter. 
Although this affects the detail structure of jet, it does not 
change the overall structure significantly. Likewise the effect 
of test section walls is not negligible, but it is not so important 
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for this experiment as long as the momentum ratio is small. 

Camera and Apparatus for Illumination 

Pictures were taken with a Polaroid camera with close-up lens. 
Its fastest shutter speed was 1/300 of a second which was a little 
too slow for taking pictures of nearly instantaneous structure of 
the jet from the side, but fast enough to see the cross-sectional 
structure . 

For illumination of the jet cross-section, an apparatus was 
made by putting a 1000 W projection lamp in the box with a slit 
in front. The slit gap was adjustable so that the gap can be made 
wide enough to give sufficient illumination yet narrow enough to 
be able to see the cross-sectional view. 

3. Experimental Procedure 

Because of the limited space in the wind tunnel test section, 
the experiment was conducted only for the case ^ = 4 =• ) • 

When side views were taken, the free stream velocity (y/p was kept 
at 5 ft /sec. The jet was illuminated from the top, and the camera 
was put about 3 feet from the side wall. The method of illumination 
and camera position for taking cross-sectional pictures are shown 
in Fig. 64. The velocity (J& was 10 ft/sec in this case. The 
slit gap and the shutter speed were determined by trial and error, 
which was made easy by the use of the Polaroid camera. The amount 
of smoke injected was about 5% of the total jet flow rate. This 
is thick enough to see the smoke clearly, yet thin enough not to 
cause change of the jet structure. 
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4. Results 

Fig. 65 shows the general structure seen from the side. A close 
look at the region near the nozzle is shown in Fig. 66. Cross- 
sectional pictures were taken at = 2, 5, 20, 70. They are 

shown in Figs. 67 through 70. Two pictures at each position are 
shown in order to see the difference in the jet structure at 
different instant. 

For the purpose of observing the vortex motion in the jet more 
clearly, a picture was taken of a laminar jet which is shown in 
Fig. 71. Since the jet velocity was too small, no data were 
taken for this case. 

5. Discussion of Results 

In the previous chapters we have seen that the jet structure 
determined from velocity, temperature or pressure distributions 
are different from each other. Therefore pictures taken here 
should be interpreted as showing the jet structure drawn from the 
concentration distribution of kerosene particles. Compared with 
the temperature distribution at the corresponding position, it can 
be seen that the concentration distribution is very similar to 
the temperature distribution. 

Fig. 65 shows the convoluted boundary between turbulence and 
non-turbulent region. In spite of the narrow space in the test 
section, the jet still has upward motion at the end of the test 
section. 

Fig. 66 gives a close look at the potential core and the 
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initial region of the jet. Before the end of the potential core, 
the jet is already deflected. The length of the potential core was 
found to be 1.7 p from the velocity measurement. 

At P/ 'p = 2 the jet cross-section is flattened and has a 
crescent shape (Fig. 67). A small amount of smoke goes into the 
separated wake. At ^/p ~ 5 the size of the jet is increased 

very much, and two distinct wakes are visible (Fig. 68). According 
to the heated jet experiment, contours of constant temperature 
at % = 20 (Fig. 46) are nearly circular; the smoke picture 

(Fig. 69) shows a distinct hollow at the bottom. The reason for 
this is that the instantaneous structure of the jet is not always 
symmetric with respect to the center-plane, as can be seen by 
comparing two pictures at each section. The symmetry plane inclines 
in a random manner, so that low temperature (low concentration) 
in the hollow region is averaged out and does not show up in the 
mean distribution. The cause of this random inclination may be 
the unbalance in the strength of two vortices. At */p ~ 70 

(Fig. 70) two vortex regions are clear, and the shape is almost 
similar to the one at ^/P = 20. 

Unlike other pictures the vortex motion shown in Fig. 71 is 
laminar. The jet is bifurcated as it goes downstream. Such 
bifurcation did not occur in the main experiment because of the 
strong interaction between two vortex wakes. 

In this flow visualization experiment, a question arises con- 
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cerning the behavior of kerosene particles which have larger size 
and weight than air. How much can they follow the streamlines of 
air? Considering horizontal motion with gravity acting downward, 
the kerosene particles have steady downward motion, the velocity 
of which can be estimated from the balance between the gravity 
force and the upward drag. Using the Stokes law of drag acting 
on a sphere in a slow uniform stream, one obtains the relation 

- 3 -ZR! (f-£) £ =■ &7Cfl) Z 

where *0 is the average radius of kerosene particles, and 

are densities of air and kerosene, respectively, and j/ is 
the steady downward velocity (the terminal velocity). Then / is 
give by 



Substituting - 10 microns and suitable values for , Jbi (of 

air) and * , the downward velocity 1/ turns out to be 3.3 x 10 2 
ft/sec, which is negligible compared with either U 0 or the average 
rotational velocity in the jet. Along the curved streamline the 
effect of the centrifugal force is important. When a kerosene 
particle moves along a curved streamline it suffers larger centrifugal 
force than air with the same velocity, so that it does not follow 
the streamline exactly. Assume that the kerosene particle deviates 
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from the streamline with the average velocity |/ / normal to the 
streamwise velocity \J . Again using the Stokes drag law, we have 

% TCR?PU' 


= E - tXM R ° V 


(37) 


where ^ is the total pressure force acting on the particle, and 
/ 

/ is the radius of curvature of the kerosene particle path. Con- 
sider a lump of air which occupies the same volume as a kerosene 
particle, and has the same velocity (J . Then we have 


% tzrIp t/ 1 ' o ' 

r r 


(38) 


where ~T is the radius of curvature of the streamline. Since V 
is expected to be small, assume /"* = Y' and £ - . Then 

Eqs. (37) and (38) give 


V _ 2. ftp / r> P j EaP. U 

u ~ ? n" r ' J r ' ? f* r . 


(39) 


In the present problem, [* is minimum at the nozzle exit (neglect- 
ing the effective source) . For simplicity assuming a parabolic 
profile of the jet centerline expressed by 


the minimum value of Y* is given by 
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» 


/ 

Therefore from Eq. (39) V/n^ is given as 


V#>M __ 4-/<_ Ro f* Uj 
U “ ? Mt>R* 


(40) 


Substituting Po~ 10 microns, - 4, (J\ ; = 5 ft/sec^ p - 1/4 in., 

and appropriate values for (of kerosene) , (of air) into 

/ 

Eq. (40), Vtr\i\X turns out to l e 



^/.2x/0 


Although this is a rough estimate, it can be concluded that 
the deviation of kerosene particles from curved streamlines of 
air are negligible. Notice that the above analysis is applicable 
to a laminar case. What this result means for a turbulent case 
is that particles are convected by the mean flow field of air. 

The fluctuating velocity causes the diffusion of particles. 



VIII. SUMMARY AND CONCLUSIONS 


Several aspects of jets in a cross flow were studied mainly 
experimentally. The problems discussed herein are: the penetra- 

tion of jets into the cross flow (jet centerline location), the 
penetration of heat (temperature centerline location) , the detail 
measurement of velocity and temperature distributions at several 
jet cross-sections, the measurement of rotational velocity distri- 
butions, study of the turbulence structure of jets, the analysis 
of the entrainment rate and finally the flow visualization 
experiment. 

The measurement of the jet centerline location in both heated 
and unheated jets shows the importance of the momentum ratio in 
the experiments. The detail measurement of velocity distributions 
reveals the marked change of the jet structure with the momentum 
ratio. Such factors as the deformation of the jet cross-section 
and the development of a vortex wake behind the jet (both are 
interrelated) contribute to produce complex patterns of the 
velocity distribution. The vortex motion grows stronger with 
increasing momentum ratio. In heated jets, this vortex motion 
is also influenced by the density ratio of two fluids. As the 
rotational velocity increases, more momentum, heat and turbulence 
energy are transported into the vortex region. The vortex 
activity attains its peak in the region where the jet centerline 
has large curvature, and decays rather slowly thereafter. The 
distributions of temperature and turbulence energy are more 
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strongly influenced by this rotational velocity field than are 
the distributions of velocity. 

The structure of turbulence in the present experiments is dif- 
ferent from that in other flows in the sense that the activity 
of large scale eddies are augmented by the vortex motion, and 
thus their energy is no longer negligible compared to that of 
energy-containing eddies. To see how much this different turbulent 
structure and also the rather complicated flow configuration 
influence the entrainment process , the entrainment rate was cal- 
culated from the data obtained here. It was found that the normal 
component of the cross flow velocity affects the entrainment rate 
in a manner nearly independent of the axial component of cross 
flow. 

Several conclusions are drawn from the present study of the 
jet in a cross flow. 

1. The jet centerline location is determined solely by the 
momentum ratio. The penetration of heat is controlled 
by both the momentum ratio and the density ratio. The 
buoyancy effect does not appear to be important within 
the density differences and in the flow region studied. 

2. The velocity distributions show that the jet structure 
is more complicated than was hitherto known. The best 
previous description of the structure was in the work 
of Jordinson (Ref. 16), in which the structure was 
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determined by total pressure distributions. The jet 
structure is strongly influenced by the momentum ratio as 
well as the density ratio for the heated jet, and this 
is expected to be true also for the cooled jet. 

3. Both temperature distributions and smoke pictures show 
that the vortex motion dominates the flow field in the far 
downstream region. 

4. The jet deflection is determined mainly by the ratio of 
the momentum flux in the direction of jet injection to 
that in the cross flow direction. The momentum flux in 
the direction of jet; injection decreases slightly because 
of the pressure force on the jet but this decrease does 
not appear to be very important to the jet deflection. 

5. The increase of the entrainment rate as compared with a 
free jet is largely due to the normal component of the 


cross flow velocity. 
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Figure 1. Coordinate System 
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Wind Tunnel Test Section 



Figure 3» Diagram of Test Apparatus 
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Figure 4. Structure of Hozzle 







Distribution of Velocity and Turbulence Intencity 
along Centerline 



Figure 6. Behavior of Free Jet from 1/4 in. Nozzle 
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Figure 7. Air Supply System 
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Figure 8. Measuring Probes 
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Figure 9. Constant-Current Hot Wire Set 





Figure 11. Diagram of Flow Direction Coordinates 
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Figure 12. Normal and Parallel Velocity Components 
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Figure 13. Location of Jet Velocity Centerline 
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Figure 14 . Correlation of Jet Velocity Centerline Location 




Figure 15* Contours of Constant Velocity 
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Figure 20. Projection of Velocity Vectors on Normal Plane 
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Figure 21* Projection of Velocity Vectors on Normal Plane 




Figure 22* Rotational Velocity Field in Vortex Region 






Distribution of Velocity Vectors 



Figure 23. Rotational Velocity Field in Vortex Region 






Figure 24 • Rotational Velocity Field in. Vortex Region 





Distribution of Velocity Vectors 



Figure 25. Rotational Velocity Field in Vortex Regio 

















Figure 29 • Mass Flux in Jet 
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Figure 33<> Static Pressure Distribution 
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Figure 34* Distribution of Turbulence Intensity 
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Figure 36* Energy Spectra of Turbulence in Jets^ 
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Figure 37 • Centerline Location of Heated and Unheated Jets 
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Figure 33. Location of Temperature Centerline 
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Figure 39. Location of Temperature Centerline ( R* 8 ) 




Figure 40. Contours of Constant Temperature 









Figure 42» Contours of Constant Temperature 







Figure 43« Contours of Constant Temperature 









Figure 45. Contours of Constant Temperature 
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Figure 46* Contours of Constant Temperature 
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Figure 49 • Temperature Distribution along Temperature Centerline 
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Figure 52. Diffusion of Temperature in Jets 




2.5 


139 



H 


iri'o- 


a 

z 


Figure 54. Comparison of Jet Centerline Location 
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Figure 55 • Correlation of Temperature Centerline Location 
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Figure 57. Potential Core Length of Jet in Cross Flow 
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Figure 59» Distribution of Entrainment Velocity 




Figure 60. Mass - Momentum Flux Diagram 
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Figure 63* Smoke Generator 
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Figure 64. Photographic Arrangement 
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Figure 65. General View of Jet in Cross Flow 



Figure 66. Close-up View 





151 



Figure 67. Jet Cross-section 
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Figure 68. Jet Cross-section 
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Figure 69. Jet Cross-section 
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Figure 70. Jet Cross-section 
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Figure 73. Z - Momentum Flux Distribution 


























APPENDIX A 


Pressure Force Acting on Jet in Cross Flow 

There are two important effects which the jet gives to the 
cross flow. One is the entrainment effect, and the other is the 
blockage effect. The former has been explained in detail in 
Chapter VI. The latter situation is similar to the case in which 
a solid object is put in a uniform flow. The flow is disturbed 
by such an object, and the object itself experiences an aerodynamic 
drag. In the present experiments, the drag in the cross direction 
acts to deflect the jet, whereas the force in the -direction 

(the direction of jet injection) reduces the jet initial momentum. 

As has been explained in Chapter VI, the jet is deflected mainly 
by the entrainment effect in the far downstream region; however 
in the region near the nozzle the pressure drag is also important. 
The decrease of the jet initial momentum can be seen in Fig. 30. 

In order to explain the effect of the pressure force analytically, 
a simplified, two-dimensional potential flow model is considered 
here. Since the cross-section is very much flattened, a two- 
dimensional model seems to be not too bad an approximation. 

Consider the potential flow around a parabolic cylinder which 
is set normal to the flow as shown in Fig. 72. Using the coordi- 
nate system in the figure, the cylinder is expressed by the equation 

-X- = -Tec 1 4r 9 ' 

or , (AD 

z = 4-az. 
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where 2a is equal to the radius of curvature at 6 =0. The 

velocity potential is given by 



The pressure distribution on the surface of the cylinder is then 

P-Pi = irfU* <*/■£&' 

(A3) 



_A_ 


where is the pressure of the undisturbed flow. It is assumed 

that the nearly same pressure is acting on a suitably defined 
surface of the three-dimensional jet. To calculate the total 
pressure force acting in the negative jZ -direction, it is neces- 
sary to consider the projection of the jet on the -plane. 
According to the flow visualization experiment, the projected area 
has nearly a rectangular shape. Thus assume that the area is 
given by 


sw = cx 


(A4) 


where the constant C represents the width of the rectangle, 
and its value is slightly different for the different momentum 
ratio. Therefore the total pressure force in the negative 
Z -direction is 
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£ = f */>-/>*£*(& 

'O * 

where is the pressure in the wake region. Generally, / “ 

Then the jet initial momentum decreases as 

Mz(x) = Mjo E 

* 


= Mj'-kfvUf 

'O 


cn+x 


where Mjo = %TCP*fC / /. 

According to the profiles of the jet centerline 

“ ~ a . 


(A5) 


Consequently 


Mjo 




U_ x \ 
8* PS 


(A6) 


Eq. (A6) is compared with the data in Fig. 73. 

As can be seen in the figure the change of the vertical momen- 
tum flux nearly follows Eq. (A6) . For £ = 4, the agreement is 
not SO good for small ^ but this is mainly due to the method of 
calculating Mz • The value of bas been obtained by inte- 

grating fj) 2 in the arbitrarily defined jet region which is a 
rather small part of the whole turbulent region for = 4 except 
in the region near the nozzle. 



APPENDIX B 


List of Symbols 

A ; 

C ; 

P 

Pt ; 

E; Et; E '2- ; 

E(f) 

Eh ; 

Q ; 

K; K / ; 

M($) ; 

Mb ; 

Mjo ; 

P Fs S 

Pe 

/*/, />* ; 

Ph ; 

£ ; 

Re ; 

SW : 

St ; 

S/2. ; 

r ; 

7c ; 


a constant 
a constant 

nozzle exit diameter 
eddy diffusivity 
entrainment coefficients 

spectrum function of turbulence kinetic energy 
Froude number 

transfer function of hot wire amplifier 
constants 

total momentum flux of jet at J* 

total momentum flux of jet at reference point 

total momentum flux of jet at nozzle exit 

pressures 

total pressure 

pressures measured by yaw meter 
Prandtl number 

? ? I 

effective velocity ratio 2 ( (J ° )2 

Reynolds number 

projected area of jet on -plane 
cross-sectional area of turbulent region 
area occupied by the contour Q =0.5 
temperature 

maximum temperature in the plane of symmetry 
at a given jet section 
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hm ; 

U ; 

Ub 


Uc 

Ue 

Urn 

; 

Ulmers ; 

Un s 

Uon ’ 

Up ' 

V, V' ; 

a ; 

fc> ; 

C L 

/ 

; 



vri 

Wt 


maximum temperature at a given jet section 
velocity 

velocity at which turbulence-non-turbulence 
interface is moving 

maximum velocity in the plane of symmetry 
at a given jet section 

entrainment velocity 

velocity component in X ^ - /*/ <2. 

maximum velocity at a given jet section 

measured value of velocity by hot wire 

velocity component - normal to jet section 

normal component of cross flow velocity 
parallel velocity components 

normal velocity components of kerosene 
particle 

a constant 

a constant 

specific heat of air at constant pressure 
mean square output voltage of wave analyzer 
frequency 

filter bandwidth of wave analyzer 
gravitational acceleration 
heat conductivity of air 
total mass flux of jet 

total mass flux of jet at reference point 
average mass flux in turbulent region 
a constant 






u,v:iir 

Uj tin 


fry ,*) 

Zo 

cX / 

P 

Y 

r ' 

<T 



S' 

0 


unit vector normal to jet section 

values of ^ (measured from temperature 
center in the plane of symmetry) at which 

temperature is equal to 1/2 ( + 7~ 0 ) 

twice kinetic energy of turbulence 


1 

•half value radius = < Sji /70 > 2 

effective radius of the region occupied by 
the contour = 0.1 

effective radius of turbulend region 

sensitivity of hot wire 

turbulence velocities 

dimensionaless velocities 

_U — Uo 

IA rr 


77 rr lSJL 

^W" Uo J n 


Cartesian coordinate system 

effective source of jet deflection 

flow direction in XZ -plane 

flow direction in -plane 

intermittency factor 

mean flow direction with respect to 
Z , -plane 

angle between (Jj, and parallel component of 

U 0 

spreading angle of mixing layer of jet in 
cross flow 

difference between £ and £ 0 

angle between yaw meter and mean flow direc- 
tion 



e 


deflection angle of jet velocity centerline 


3 

ft 

v 

X 

f 

f' 

v 

(i () 


dimensionless temperature 

ZzJL 

dynamic viscosity 
kinematic viscosity 
location of jet virtual origin 
density 

density fluctuation 
similarity variable 

energy/cps of wave analyzer input voltage 
jet defined coordinate system 


( ) ; 
subscript ; 

subscript ; 


time-mean value 

quantity of undisturbed cross flow 
quantity of initial jet 
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